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Me  analyze  the  theoretical  basis  for  th£fmb:M^and  Variable  Frequency  Mag¬ 
nitude  (VFM)  earthquake/explosion  discriminants  in  light  of  recent  results  from 
nunerlcal  modeling  of  earthquake  and  explosion  dynamics.  Me  examine  the  effects 
of  source  mechanism  (double-couple  versus  point  dilatation),  source  depth, 
source  spectra,  and  source  region  elastic  structure.  We  use  six  three-dimension¬ 
al  finite  difference  slnulatlons  of  earthquakes  for  the  earthquake  source  models, 
and  four  empirical  (Mueller -Murphy)  and  four  numerical  models  for  explosions  — 
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ABSTRACT:  (continued) 

in  different  materials.  We  estimate  that  source  mechanism  (i.e.,  quadrupole 
versus  monopole)  contributes  about  .35  magnitude  on  the  average  to  the  sep¬ 
aration  between  earthquakes  and  explosions  on  a  mb:Ms  diagram.  >  Differences 
in  source  region  elastic  properties  tend  to  increase  the  separation,  contrib¬ 
uting  a  factor  of  approximately  /pa3  (p  =  density,  a  =  P-wave  speed)  to  the 
ratio  of  surface  wave  amplitude  to  body  wave  amplitude.  For  events  with 
Ms  <  4.0,  this  effect  may  be  more  important  than  source  spectral  differences, 
especially  for  explosions  in  low  velocity  materials. 


VFM  discrimination  depends  primarily  on  spectral  differences  between 
earthquakes  and  explosions.  Theoretical  separation  of  the  earthquake  and 
explosion  populations  is  very  clear  for  explosions  in  granite  at  frequencies 
\  Hz  or  greater  with  m^  >  4.0.  Smaller  events  may  be  discriminated  using 
higher  frequencies.  Separation  is  not  as  good  for  explosions  in  low  veloc¬ 
ity  materials  such  as  tuff,  especially  at  low  yields. 

For  small  events,  the  mb:M$  and  VFM  discrimination  methods  work  (or 
fail)  for  quite  different  reasons.  Interference  with  pP  causes  earthquake 
and  explosion  populations  to  converge  on  the  mb:Ms  plane  at  mh  4.5  using  an 
empirical  depth/yield  relation  of  d  =  122W1/3  meters.  This  interference  in¬ 
creases  the  separation  using  VFM,  allowing  the  method  to  work  even  for  very 
low  yields.  More  deeply  buried  small  explosions  may  be  discriminated  by 
mb;Ms ,  but  not  by  VFM.  Simultaneous  use  of  both  methods  should  allow  im¬ 
proved  discrimination. 
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I.  INTRODUCTION 


A  number  of  methods  have  been  proposed  for  discriminating 
earthquakes  from  underground  nuclear  explosions  using 
characteristics  of  the  seismic  signal  observed  at  teleseismic  and 
regional  distances.  The  mb:M$  and  Variable  Frequency  Magnitude 
(VFM)  discriminants,  in  particular,  have  received  extensive  study, 
and  each  has  a  wel 1-estaol ished  empirical  basis.  Theoretical 
understanding  of  these  discriminants,  however,  is  quite  incomplete. 
In  view  of  the  prominent  role  envisioned  for  seismic  methods  in 
detecting  ana  identifying  underground  nuclear  explosions,  it  is 
important  that  these  techniques  be  quantitatively  understood  on  the 
oasis  of  eartnquake  and  explosion  physics.  Only  with  a  quantitative 
model  describing  why  the  discriminants  work,  can  we  confidently 
predict  the  circumstances  under  which  they  may  fail. 

The  difference  between  body  wave  magnitude,  mb,  and  surface 
wave  magnitude,  M$,  is  usually  larger  for  explosions  than  for 
earthquakes.  The  separation  between  the  earthquake  and  explosion 
populations  on  the  mb:Ms  plane  has  been  documented  in  many 
observational  studies.  Marshall  and  Basham  (1972),  for  example, 
show  m^  separation  for  Eurasian  events  with  a  given  M ^  ranging 

from  aDcut  0.5  to  1.5  mb  units  with  somewhat  better  separation  for 
1 arger-magn i tude  events  than  for  small  ones.  Another  proposed 
discriminant,  the  Variable  Frequency  Magnitude  (VFM)  method 
(Arcnambeau,  Harkrider  and  Helmberger,  1974)  is  based  on  the 

spectral  content  of  the  short-period  P-wave.  In  this  method,  the 
body-wave  magnitude  is  measured  from  narrow-band-filtered 
seismograms  at  two  different  frequencies,  f1  and  f2-  A  plot  of 

mb(fi)  versus  mQ(f2)  also  produces,  in  many  circumstances,  a 
clear  separation  of  the  eartnquake  and  explosion  populations 
(Savino,  Archambeau  and  Masso,  1980). 

Altnough  the  empirical  results  are  well-known,  theoretical 
understanding  of  the  discriminants  has  been  impeded  by  lack  of 

adequate  models  of  the  earthquake  source.  Kinematic,  semiempirical 


models  of  eartnquakes  such  as  the  Haskell  (1964,  1966),  Brune 

(1970),  and  Aki  (1967,  1972)  models  have  been  available,  but  these 
models  contain  free  parameters  such  as  the  rise  time  or  corner 

frequency  which  are  difficult  to  constrain  and  crucial  to  the 
analysis  of  discrimination  methods.  Aki,  Bouchon  and  Reasenberg 
(1974)  apply  the  earthquake  spectral  theory  of  Aki  to  argue  that 

source  spectrum  (or  equivalently,  source  time  history)  alone  can 

account  for  about  1  magnitude  unit  separation  in  mb  for  events 

with  M  as  low  as  2.  On  the  other  hand,  Douglas,  Hudson  and 
Kemohavi  (1971),  using  point  sources  in  a  half-space,  argue  that 

differences  in  source  mechanism  (double-couple  versus  point 
dilatation)  alone  are  sufficient  to  give  reasonable  discrimination 
between  explosions  and  earthquakes  except  for  a  few  preferred 
orientations  of  tne  double-couple.  Douglas,  Hudson  and  Blarney 

(1972),  using  synthetic  seismograms  for  layered  media,  and  using 
point  and  extended  earthquake  sources,  come  to  the  conclusion  that 
source  mechanism  is  most  important  to  discrimination.  They  also 

find  that  explosions  in  soft  materials  have  a  greater  separation 

between  m.  and  M  than  explosions  in  harder  materials. 

D  S 

Archambeau,  Harkrider  and  Helmberger  (1974)  first  developed 
tne  physical  basis  for  the  VFM  discriminant  and  explained  the 
differences  between  m^  and  of  earthquake;  and  explosions  on 

the  basis  of  characteristic  spectral  differences.  They  also  looked 
at  applications  of  the  VFM  method  to  the  discrimination  of  multiple 
explosions.  The  earthquake  model  used  in  that  study  was  a  growing, 
propagating,  transparent  spherical  cavity  (Archambeau,  1968).  The 
explosion  model  used  was  an  early  version  of  the  finite  difference 
models  used  here.  In  this  study,  we  reduce  uncertainties  about  the 
source  spectra  by  using  a  more  appropriate  failure  geometry  with 
more  realistic  boundary  conditions  for  our  earthquake  models,  more 
refined  numerical  explosion  models,  and  an  empirical  explosion  model 
which  is  constrained  by  near  field  observations. 

Von  Seggern  and  Rivers  (1979)  evaluated  the  reliability  of  a 
number  of  earthquake/explosion  discriminants,  using  methods  similar 
to  those  used  in  tnis  report,  and  applied  them  to  events  in  the 
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western  United  States.  In  general,  their  conclusions  agree  with 
ours;  however,  there  are  three  important  differences.  First,  our 
improved  source  models  allow  increased  confidence  in  our  estimates 
of  spectral  differences  between  earthquakes  and  explosions.  Second, 
in  their  estimation  of  P-wave  amplitudes,  they  use  a  whole  space 
solution,  neglecting  the  transmission  coefficient  (Equation  (5.2)) 
which  shoulu  be  included  if  the  earthquake  and  explosion  are  in 
different  materials.  This  transmission  coefficient  reduces  the 
separation  between  earthquakes  and  explosions  on  the  mb:Ms 
plane.  Finally,  their  estimate  of  source  mechanism  effects  was 
basea  on  root  mean  square  averages  of  the  radiation  pattern  while  we 
computed  distribution  functions  by  directly  calculating  amplitudes 
for  a  large  number  of  orientations. 

t  It  is  important  to  understand  the  relative  contribution  of 

source  mechanism  and  source  spectrum  to  and  of  earthquakes 
since  the  two  phenomena  have  very  different  implications  for  the 
applicaoility  of  tne  discriminant  at  low  magnitude.  The  effect  of 
t  source  mechanism  is  independent  of  source  size  whereas  the  effect  of 

spectral  shape  inevitably  vanishes  below  some  threshold  size. 
Furthermore,  if  the  observed  mb-Ms  earthquakes  is 

predominantly  an  effect  of  source  spectrum  rather  than  source 
mecnanism,  then  the  discriminant  will  be  more  likely  to  fail  for 
high  stress-drop  earthquakes. 

Conclusions  about  the  influence  of  earthquake  complexity  on 
mb:Ms  are  also  highly  dependent  on  the  source  model.  31andford 
t  (1975)  suggested  tnat  increasing  the  complexity  of  an  earthquake,  by 

subdividing  into  many  small  earthquakes,  causes  an  ever  increasing 
earthquake  corner  frequency.  This  would  imply  that  increasing  the 
complexity  makes  the  earthquake  more  explosion-like. 

The  purpose  of  tnis  study  is  to  refine  the  theoretical  basis 
of  the  n>b:Ms  and  VFM  methods.  Clearly,  it  is  important  for  this 
purpose  to  use  an  earthquake  source  model  which  is  analytically 
sound  and  relatively  free  from  non-physical  assumptions.  We  will 
use  theoretical  source  spectra  (and  time  histories)  obtained  from 
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three-dimensional  numerical  simulations  of  earthquakes  by  Day  (1979, 
1982a,  1982b).  A  range  of  sources  is  considered,  including  uniform 
stress-drop  events,  isolated  stress  concentration  events,  and 
multiple-stress-concentration  events.  We  use  two  types  of  explosion 
sources  for  comparison:  the  one-dimensional  (spherically  symmetric) 
numerical  simulations  of  Cnerry,  Rimer  and  Wray  (1975),  and  the 
empirical  explosion  source  model  of  Mueller  and  Murphy  (1971). 

We  first  give  separate  consideration  to  each  of  four  factors 
likely  to  contribute  to  differences  between  earthquake  and  explosion 
seismograms.  These  are  source  mechanism  (i.e.,  double-couple  versus 
center-of-dilatation),  source  depth,  source-region  elastic 
structure,  and  source  spectrum  (or  time-history).  After  analyzing 
tnese  factors  separately,  we  make  time-domain  mb  estimates  using 
synthetic  seismograms  which  incorporate,  simultaneously,  effects  of 
free-surface  reflection,  crustal  structure,  source  mechanism,  and 
source  finiteness.  We  restrict  consideration  to  earthquakes  and 
explosions  occurring  in  the  same  general  area  so  that  travel  path 
and  attenuation  effects  can  be  considered  the  same  for  all  sources. 
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II.  SOURCE  MODELS 


2.1  EARTHQUAKE  MODELS 

Dynamic  earthquake  models  in  the  sense  used  here  are  those 
which  specify  the  tectonic  prestress  and  the  resisting  (frictional) 
stresses  in  the  source  region,  and  derive  the  resulting  displacement 
discontinuity  (slip)  via  the  equations  of  motion.  Such  earthquake 
models  generally  lead  to  nonlinear,  mixed  boundary  value  problems 
requiring  numerical  methods  for  their  solution.  Only  recently  have 
three-dimensional  solutions  for  finite  faults  become  availade.  The 
dynamic  models  for  which  numerical  solutions  have  been  obtained  fall 
into  two  catagories:  f ixed-rupture-velocity  models,  in  which 

rupture  velocity  is  a  specified  constant,  and  spontaneous  rupture 
models,  in  which  rupture  velocity  is  determined  by  some  failure 
criterion,  then  derived  as  part  of  the  numerical  solution. 

Six  of  tne  aynamic  earthquake  simulations  performed  by  Day 
(1979,  1982a, b)  provide  the  basis  for  the  earthquake  source 

functions  usea  in  this  study.  Three  of  these  are  for  fixed-rupture- 
velocity  models  and  three  are  for  spontaneous  ruDture  models.  The 
former  are  uniform  stress-drop  events,  while  the  latter  incorporate 
various  forms  of  non-uniform  stress-drop.  The  physical  basis  of 
tnese  simulations,  as  well  as  a  detailed  discussion  of  the  resultant 
fault  slip  histories,  is  given  in  the  above  references.  Here  we 
only  outline  the  main  characteristics  of  the  computed  rupture 
histories  and  far-field  seismic  radiation. 

The  final  perimeters  of  the  six  fault  surfaces  are  sketched  in 
Figure  1.  Models  3,  4,  and  5  are  fixed-rupture-velocity  cases  in 
which  the  rupture  grows  circularly  outward  from  the  center  of  the 
fault  at  a  velocity  of  90  percent  of  the  shear  wave  speed  s,  and 
stops  abruptly  at  the  prescribed  perimeter.  Models  1,  2,  and  6  are 
spontaneous  rupture  cases.  For  these  models,  the  distribution  of 
prestress  on  tne  plane  of  faulting  is  as  shown  in  the  figure.  The 
dynamic  stress-drop,  which  is  the  difference  between  the  shear 
prestress  and  the  sliding  frictional  stress  on  the  fault,  is 
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Figure  1.  Sketches  of  six  earthquake  source  models  used  in  this  study. 

The  final  half-length  of  the  fault  is  denoted  by  a,  the  maximum 
dynamic  stress-drop  by  Ax.  Models  3,  4,  and  5  are  fixed  rupture 
velocity  models  and  1,  2,  and  6  are  spontaneous  rupture  models. 
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prescribed  to  vary  over  the  plane  of  faulting  as  shown  in  the 
figure,  witn  at  representing  its  maximum  value.  Average  dynamic 
stress-drop  in  these  models  turns  out  to  be  nearly  the  same  as 
average  static  stress-drop,  which  is  an  often-quoted  seismic  source 
parameter. 

Models  i  and  2  represent  isolated  stress  concentrations,  while 
Model  6  represents  an  assemblage  of  five  discrete  stress 
concentrations.  The  final  fault  perimeters  for  Models  1  and  2  were 
not  prescribed  in  the  calculations;  rather,  the  rupture  front 
aecelerated  and  stopped  at  the  position  shown  in  response  to  the 
prescribed  stress  inhomogeneity.  For  Model  6,  rupture  growth  was 
stopped  by  a  prescribed  jump  in  shear  strength  at  the  rectangular 
perimeter  shown  in  the  figure. 

Rupture  velocity  turns  out  to  be  highly  variable  for  each  of 
the  Models  1,  2,  and  6,  as  discussed  by  Day  (1979,  1982b).  For 
example,  in  eacn  case,  the  average  rupture  velocity  (final 
half-length,  a,  divided  by  rupture  arrival  time  at  the  ends  of  the 
fault),  whicn  is  given  in  the  figure,  is  less  than  a;  but  in  each 
case,  local  rupture  velocity  somewhere  on  the  fault  exceeds  s. 

These  six  earthquake  models  incorporate  a  fairly  wide  range  of 
pnenomena:  fault  length-to-width  ratios  between  1  and  16/3,  uniform 
versus  localized  prestress,  and  constant  versus  variable  rupture 
velocity.  There  is  consideraDle  observational  evidence  that 
variaole  stress-orop  effects  are  prevalent  in  shallow  earthquakes 
(for  example.  Hanks,  1974;  Hartzell  and  Brune,  1979;  and  Sache, 
Lamoert  and  Barker,  1980),  ana  the  variety  of  source  moaels  used  in 
tnis  study  will  permit  us  to  make  at  least  qualitative  assessments 
of  their  importance  for  seismic  discrimination.  This  is  the  concern 
addressed  Dy  Blandford  (1975),  using  an  ad  hoc  model  for  the  source 
spectrum.  Our  results,  based  on  the  above  dynamic  models,  do  not 
agree  with  his  conclusions.  In  fact,  for  the  models  considered 
here,  source  complexity  tends  to  reduce  the  source  spectrum  near  the 
corner  although  this  reduction  may  be  accompanied  by  spectral 
enhancement  well  above  the  corner. 


Figure  2  shows  the  far-field  P-  and  S-wave  displacement 
spectra  and  time  histories  radiated  by  these  six  fault  models.  For 
a  given  fault  orientation,  the  seismic  radiation  depends  on  both 

azimuth  and  takeoff  angle,  and  the  results  in  Figure  2  are,  in  each 
case,  for  a  45  degree  dip-thrust  orientation  observed  at  azimuths  of 
+•  90  degrees  from  the  fault  strike.  Takeoff  angles  are  appropriate 
to  excite  teleseismic  P-waves  at  a  range  of  about  35  degrees  (»  * 
27*  for  P-waves,  s  =  15*  for  S-waves). 

Also  shown  in  Figure  2  are  the  RMS-averaged  (over  e  and  t) 

spectra  for  P-  and  S-waves.  The  RMS  spectra  are  much  simpler  than 

the  individual  spectra  since  the  focal-spnere  averaging  has  removed 

spectral  holes  present  for  individual  rays.  They  indicate  that 

differences  in  average  spectral  shape  between  P-  and  S-waves  are 

surprisingly  small,  particularly  for  the  spontaneous  rupture  cases. 

Their  smoothness  makes  the  RMS  spectra  particularly  useful  for 

making  spectral  estimates  of  body-wave  magnitude  as  is  done  in  a 

subsequent  section.  On  the  other  hand,  actual  body-wave  magnitudes 

are  time-domain  measurements  made  from  a  short  segment  of  the 

wavetrain,  and  pnase  effects  may  be  quite  important,  particularly 

for  complex  waveforms  such  as  those  raaiated  by  Model  6.  Therefore, 

we  will  also  use  the  individual  (amplitude  and  phase)  spectra,  in 

comoination  with  synthetic  seismogram  methods,  to  maxe  time-domain 

m.  estimates  for  the  sources, 
b 

The  spectra  and  time  histories  in  Figure  2  are  snown  in 

non-dimensional  form.  When  we  construct  spectral  and  time-domain 
magnitude  estimates,  it  will  be  necessary  to  scale  each  source 
function  to  various  values  of  the  fault  half-length,  a,  and  maximum 
stress-drop,  at.  The  assumptions  implicit  in  such  scaling  are 
described  by  Day  (1982b).  In  short,  we  must  assume  that  the 

following  dimensionless  quantities  are  invariant  with  source  size 
and  stress-drop:  1)  Poisson's  ratio  of  the  source  medium,  2)  the 
ratio  of  fault  shear  strengtn  to  dynamic  stress-drop,  and  3)  the 
ratio  of  total  fracture  energy  to  total  strain  energy  release. 
Under  these  assumptions,  the  rupture  history  remains  constant  as  the 
dimension  of  the  fault  is  changed. 
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As  indicated  in  Figure  2,  spectral  ordinate  varies  directly 
witn  stress-drop  and  fault  length  cubed,  and  inversely  with 
wavespeed  of  the  source  medium  cubed,  while  frequency  scales 
inversely  with  length  and  directly  with  the  wavespeed  of  the  source 
medium.  In  subsequent  analyses,  we  will  set  the  P-  and  S-wave 
speeds  of  the  earthquake  source  region  to  6  km/sec  and  3.5  km/sec, 

3 

and  tne  density  to  2.7  gm/cm  ,  values  appropriate  for  granite  at 
mid-crustal  depths.  Unless  otherwise  noted,  we  will  set  a  to  100 
bars  (iO  MPa).  This  value  implies  average  stress  drops  ranging  from 
43  oars  to  100  bars  as  shown  in  Figure  1.  With  the  velocities  and 
stress  drop  held  constant,  a  simple  cube  root  scaling  applies: 
multiplying  the  spectral  amplitude  by  a  factor  A  multiplies  the 
corresponding  frequency  by  A-^. 

2.2  EXPLOSION  MODELS 

We  consider  two  types  of  explosion  source  models.  The  first 
are  one-dimensional  (spnerically  symmetric)  finite  difference 
simulations  of  explosions  buried  in  different  geologic  materials. 
The  method  of  computation  is  described  by  Cherry,  Rimer  and  Wray 
(1975).  The  calculations  used  here  were  performed  more  recently  by 
N.  Rimer  (unpublished  worn  performed  at  Systems,  Science  and 
Software,  1978-1979)  for  explosions  in  dry  granite,  saturated 
granite,  salt  and  saturated  tuff.  The  calculation  for  salt  was 
intended  to  model  the  explosion  SALMON  (yield  5.3  KT).  The 
saturated  granite  calculation  is  a  model  for  the  explosion 
PILEDRIVER  (yield  60  KT) .  The  tuff  mod.il  is  a  typical  Nevada  Test 
Site  material.  The  dry  granite  calculaUo.i  is  a  model  for  French 
explosions  in  tne  Hoggar  region  of  the  jahara.  This  calculation 
leads  to  a  very  low  amplitude  reduced  velocity  potential  compared  to 
explosions  in  the  other  materials  at  the  same  yield. 

We  use  the  reduced  velocity  potential  y (w)  derived  from  these 
calculations  to  make  spectral  magnitude  estimates  and  time-domain 
synthetic  seismograms.  The  reduced  velocity  potential  may  be  cube 
root  scaled  to  arbitrary  yield:  the  spectral  amplitude  scales 


directly  with  yield  while  the  frequency  varies  as  the  inverse  cube 
root  of  the  yield.  The  explosion  elastic  radius  varies  directly 
with  the  cube  root  of  the  yield.  Although  this  scaling  may  be  done 
for  any  yield,  the  numerical  calculations  were  performed  for  a  fixed 
overburden  pressure.  In  general,  larger  explosions  will  be  more 
deeply  buried  than  smaller  explosions;  so  the  scaling  for  'i'  may  not 
oe  applicable  to  explosions  at  all  yields. 

For  the  second  explosion  source  model,  we  use  the  empirical 
source  model  of  Mueller  and  Murphy  (1971).  Details  of  the  model  are 
given  in  Appendix  A.  This  model  includes  a  scaling  with  depth  as 
well  as  a  scaling  with  yield  (W).  The  main  result  of  the  depth 
dependence  is  that  tne  explosion  elastic  radius  varies  approximately 
as  W1'5  rather  than  because  of  the  increasing  overburden 

pressure.  We  have  used  Mueller-Murphy  reduced  velocity  potentials 
for  explosions  in  salt,  granite,  shale  and  tuff/rhyol ite. 

Table  1  gives  the  material  properties  of  the  source  media  for 
the  eight  explosion  sources  (a  =  congressional  wave  velocity,  s  = 
snear  wave  velocity,  p  =  density).  The  long  period  limit  of  the 
reduced  velocity  potential,  'i'oo,  and  the  moment  MQ  (defined  by  MQ 
=  4woa2 'i'a,)  are  given  for  each  source  scaled  to  100  KT.  For  the 
Mueller-Murphy  sources,  depth  and  overburden  pressure  are  a  function 
of  yield;  for  tne  finite  difference  calculations  they  are  fixed. 
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III.  EFFECT  OF  SOURCE  MECHANISM 


As  a  first  step  in  comparing  earthquakes  and  explosions  as 
generators  of  body  (P)  and  surface  (Rayleigh)  waves,  we  estimate 
differences  in  body-wave  magnitude  (m^)  and  surface-wave  magnitude 
(M$)  attributable  solely  to  the  difference  between  double-couple 
and  center-of-dilatation  point  sources.  We  neglect,  in  this 
comparison,  any  characteristic  differences  in  source  spectrum  and 
source-region  materi al  properties. 


In  a  uniform,  infinite  space,  we  can  write  the  spectral 
amplitude  of  the  far-field  P-wave  radiated  by  an  earthquake, 


A^  (id, 9,0) 


mn(<*>,  9.  0) 

-3 - 3 -  R„(®.  *) 

4irpa  r 


(3.1) 


where  R  (o,  0)  is  the  double-couple  radiation  pattern  factor  (with 
maximum  absolute  value  of  1),  m  (u,  9,  0)  is  the  earthquake 
source-spectrum  or  moment  function,  p  is  the  density,  a  is  the 
P-wave  speed,  and  r,  9,  and  0  are  the  spherical  coordinates  of  the 
receiver  relative  to  the  hypocenter.  The  source-spectrum  mq(u>,  9, 
0)  contains  all  the  effects  of  spatial  and  temporal  finiteness  of 
the  earthquaxe  (Kostrov,  1968).  The  long-period  limit  of  the  source 
spectrum,  which  is  indepenaent  of  9  and  0,  is  the  seismic  moment  of 
tne  earthquake,  denoted  M^. 

Similarly,  the  far-field  P-wave  spectral  amplitude  of  an 
underground  explosion,  a£,  is  given  by  (assuming  spherical 
symmetry  of  tne  source) 


mx(«) 

4irpa  r 


(3.2) 


where  mx(w)  is  the  explosion  source  spectrum.  The  long-period 

limit  of  the  source  spectrum  is  the  seismic  moment  of  the  explosion, 

Mx.  The  source  spectrum  m  (id)  is  proportional  to  the  reduced 
0  x  # 

velocity  potential  of  the  explosion,  'J'(id): 


4irpa^f(u) 


(3.3) 
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t 


» 


t 


mx(u)  * 


We  now  define  Ap  and  A*  to  be  the  Rayleigh  wave 
spectral  amplitudes  of  the  earthquake  and  explosion,  respectively. 
We  define  the  "relative  amplitude  factor,"  X,  to  be 


X 


AX 

AX 


(3.4) 


which  gives  the  generation  of  Rayleigh  waves  by  an  earthquake 

relative  to  the  generation  of  Rayleigh  waves  by  an  explosion  of  the 
same  oody-wave  magnitude.  Its  logarithm,  whicn  we  call  the 
"relative  magnitude  factor,"  is 

logio'X  =  (Mj-Mj)  -  ,  (3.5) 

which  gives  tne  difference  (Mp  -  M*)  for  an 

earthquake  and  explosion  with  the  same  mb,  or  the  mb  difference 

(m*  "  mb)  for  an  event  pair  with  the  samfe  V  In  the 

following  analysis,  we  will  be  considering  distributions  of  the 

relative  magnituae  factor  for  an  ensemble  of  events.  The  mean  value 
of  log^gX,  tnen,  will  represent  the  mean  distance  between  the 
earthquake  and  explosion  populations  on  the  mb:M$  plane. 

We  may  obtain  estimates  for  the  relative  amplitude  factor 
using  expressions  given  by  Douglas,  Hudson  and  KemDhavi  (1971).  The 
earthquake  is  represented  by  a  point  double-couple  (dislocation) 
with  fault  normal  n  and  unit  slip  vector  s,  located  at  depth  h,  and 
the  explosion  by  a  center-of-di latation  located  at  depth  h*.  If  the 
earthquake  and  explosion  have  the  same  source  spectral  shape  and  the 
same  seismic  moment  (mx(w)  *mb(iu)),  Douglas  et  al .  give  the 
following  amplitude  ratio  for  Rayleigh  waves: 
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Aj  2  <-Hn  (h*  -  h) 

-i  ■  \  «p  -5jl_ - 

AS  =R  * 


"o  31- 


2q1ina(l  -  E) 


*"2  1-1-^7  E 


where 

3  Rayleigh  wave  angular  frequency 

cR  *  Rayleigh  wave  speed 

a, 3  »  P  and  S  wave  speeds 

M  e2 1  2  \  1/  2 
na  =  "  CR/a  ' 

n  r2/82l1/2 

n3  "  t1  -  <VS  } 

9  *  Takeoff  angle 

4>  3  Azimuth  of  recording  instrument 

qo  a  s3n3 

ql  *  ( sln3  *  s3ni)cos*  +  (s2n3  +  s3n2^  sine 
q2  *  ^slnl  “  s2n2 )cos2^  +  (sin2  +  s2nl)  sin2^ 

uiRh 

E  «  e*P  ^  '  V 

L 

For  the  following  analysis,  we  take  a  *  3,  which  requires  Cp 


From  equations  3.1  and  3.2,  the  P-wave  amplitude  ratio, 
A?/A?  for  events  of  the  same  moment  ana  spectral  shape, 
is  simply  equal  to  the  modulus  of  the  radiation  pattern  factor, 
|r  |.  This  is  given  by  Douglas  et_  al .  in  the  same  notation  as 

Equation  (3.6): 


SCUBED 


(3.7) 


Ah  i  ?  2  i 

-j  =  |  qQ(2  -  3  sin  a)  +  q^.  sin2a  +  q2  sin  a  j  . 

Ab 

These  expressions  apply  to  a  uniform  half- space,  and  will  be  the 
basis  for  the  analysis  which  follows  in  this  section.  Equation  3.7 
neglects  effects  of  surf ace-reflected  waves.  Modifications 
introduced  by  reflected  phases  and  near-source  material  properties 
will  be  reviewed  in  later  sections. 

Douglas  et_  al^  used  Equations  3.6  and  3.7  to  estimate  the 
effects  of  the  source  mechanism  on  the  relative  magnitude  factor. 
Some  of  the  conclusions  of  the  Douglas  et_  ak.  paper,  however,  depend 
rather  strongly  on  the  choice  of  a  very  small  takeoff  angle  of  10 
degrees  and  on  the  limited  number  of  source  orientations 
considered.  This  takeoff  angle  is  appropriate  for  a  source-layer 
P-velocity  of  3  km/second  for  an  observation  point  located  65 
degrees  away.  This  assumption  led  to  the  conclusion  that  a  shallow, 
vertical  dip-slip  fault  is  a  more  efficient  generator  of  surface 
waves  compared  to  body  waves  than  a  45  degree  dip-slip  fault.  This 
result  arises  because  the  assumed  takeoff  angle  is  close  to  a  =  0, 
whicn  is  a  P-wave  node  for  the  vertical  dip-slip  fault  and  an 
anti  node  for  the  45  degree  dip-slip.  A  more  appropriate  crustal 

compressional  velocity  for  all  but  the  shallowest  of  earthquakes  is 
about  6  km/second.  This  requires  takeoff  angles  of  between  29 
degrees  and  17  degrees  for  an  observation  distance  range  of  30 
degrees  to  80  degrees  (Richter,  1958). 

In  TaDle  2,  we  have  listed  the  relative  amplitude  factors  for 
vertical  dip-slip,  45  degree  dip-slip,  and  strike-slip 

double-couples  using  an  explosion  depth  of  one  kilometer.  We  have 
used  a*6  km/sec,  a  period  of  20  seconds,  and  azimuths  of  90 
degrees  for  the  dip-slip  events  and  45  degrees  for  the  strike-slip 

events  (all  maxima  of  the  Rayleigh  wave  radiation  pattern).  The 

.  2 

numbers  on  the  table  must  be  divided  by  sin2a,  cos2a,  and  sin  a 
respectively  to  include  the  effect  of  the  body  wave  radiation 
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TABLE  2 
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RELATIVE  AMPLITUDE  FACTORS  FOR  POINT  DISLOCATIONS 


Earthquake 
Depth  (km) 

Vertical 

Dip  Slip 

45  * 

Dip  Slip 

Strike  Slip 

1 

.27 

1.80 

1.40 

5 

.79 

.70 

.70 

10 

1.60 

.15 

.16 

15 

1.80 

.60 

.14 

20 

1.80 

.80 

.30 

25 

1.70 

.88 

.37 

30 

1.50 

.86 

.38 

l 
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pattern.  The  tnree  sources  behave  very  differently  with  changes  in 
depth.  For  very  shallow  (<  5  km)  sources,  the  45  degree  dip-slip 
and  strike-slip  events  are  much  more  efficient  generators  of  surface 
waves  than  the  vertical  dip-slip  events.  On  the  other  hand,  for 

depths  greater  than  5  km,  the  vertical  dip-slip  becomes  a  more 

efficient  generator  of  surface  waves.  The  45  degree  dip-slip  and 
strine-slip  pass  through  a  node  and  are  very  inefficient  surface 
wave  generators  at  about  10  km  depth. 

In  order  to  avoid  biases  associated  with  specific  source 

orientations,  we  examine  the  frequency  distribution  of  the  relative 
amplitude  factor  for  each  of  several  ensembles  of  earthquaxes.  In 
each  ensemble,  it  is  assumed  that  the  distribution  of  fault 

orientations  is  uniform  in  the  sense  that  observations  at  a  given 
takeoff  angle  and  azimuth  sample  the  focal  sphere  uniformly.  Figure 
3  shows  the  resulting  distribution  ^unctions  of  the  relative 
magnitude  factor,  calculated  from  Equations  3.6  and  3.7,  for  each  of 
several  earthquake  source  depths  and  takeoff  angles.  The  explosion 
source  depth  is  1  km  in  each  case.  Recall  that  the  horizontal  axis 
represents  the  estimated  separation  in  the  :MS  plane  between 
an  explosion  and  an  earthquake  of  the  same  mb  (or  the  same  Ms) . 
The  vertical  axis  can  be  thought  of  as  the  relative  number  of 
observations  of  event  pairs  with  a  given  separation  when  effects  of 
source  spectrum,  source  region  material  properties,  and  reflected 
phases  are  neglected. 

The  distributions  in  Figure  3  are  surprisingly  smooth. 
Although  the  excitation  of  Rayleigh  waves  relative  to  P-waves  for  a 
particular  fault  orientation  is  quite  sensitive  to  source  depth  and 
takeoff  angle,  the  distribution  of  tht  relative  excitation  over  a 
range  of  fault  orientations  is  not  very  sensitive  to  source  depth  or 
takeoff  angle.  The  average  separation  of  the  populations  due  to 
source  mechanism  is  only  about  one  third  of  a  magnitude,  and  about 
one  third  of  the  observations  have  a  relative  magnitude  factor  less 


than  zero. 
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10°,  20°  and  30“  respectively. 


The  values  of  mb  and  M$  used  in  discrimination  studies  are 

usually  averages  of  mb  and  Ms  over  a  range  of  observation 

azimuths.  If  we  adopt  this  approach,  and  average  the  theoretical 

Rayleigh  wave  and  body  wave  magnitudes  separately  instead  of  looking 

at  individual  stations,  then  we  get  a  more  realistic  estimate  for 

mb  and  M$.  We  can  also  consider  a  range  of  depths  and  takeoff 

angles.  Figure  4  shows  the  distribution  function  for  event  pairs 

with  M  and  mu  averaged  over  azimuth  for  earthquakes  uniformly 
s  b 

distributed  in  source  depth  (l  to  25  km)  and  takeoff  angle  (10  to  30 
degrees)  as  well  as  in  source  orientation.  Points  very  close  to 
nodes  (less  than  5  percent  of  the  maximum  amplitude)  were  not  used 
in  the  average  of  magnitudes  over  azimuth.  The  result  is  a  narrower 
distribution  function  sharply  peaked  at  about  .3  -  .4  magnitude. 
The  average  over  azimuth  removes  tne  extremes  which  may  be  seen  in 
individual  observations  resulting  from  observation  points  close  to 
nodes  in  the  body  or  surface  wave  "adiation  patterns. 

We  can  see  the  individual  contributions  of  differences 
and  m^  differences  by  computing  and  m^  separately  using 

Equations  3.6  and  3.7  and  plotting  a  two-dimensional  distribution. 
The  contour  plot  in  Figure  5  shows  the  frequency  of  events  as  a 
function  of  M$  difference  and  mb  difference  for  earthquakes  and 

explosions  with  the  same  moment.  Again,  these  distributions  are  for 
a  range  of  depths  and  takeoff  angles  and  are  averaged  over  azimuth. 
For  all  event  pairs,  mb  of  the  explosion  exceeds  mb  of  the 

earthquake.  This  is  simply  a  result  of  the  radiation  pattern 
difference  mentioned  before.  A  more  surprising  result  is  that 
earthquakes  have  slightly  smaller  values  than  explosions  with 
the  same  moment.  The  distributions  are  strongly  concentrated  about 
the  point  *  -.1»  mb-mb  m  Fifty 

percent  of  the  event  pairs  are  located  in  the  interval  -.3  <_ 

_M*  <  0.,  -.6  <  mj-mj  <  -.3.  About  20  percent  of 

event  pairs  have  ^s-,^s""^mb_n'b^  —  an<^  cannot 

oe  discriminated  by  source  mechanism  alone. 
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Magnitude  Difference 


Figure  4.  Frequency  distribution  of  the  difference  between  earthquake  and 

explosion  magnitudes  due  to  source  mechanism  (double-couple  versus 
dilatation)  alone.  The  horizontal  axis  is  the  function  M*  (earth¬ 
quake)  -  Ms  (explosion)  -  (mb(earthquakes)  -  mb(explosionl).  For 
two  earthquakes  with  the  same  Ms,  it  represents  the  amount  by  which 
mb(explosion)  exceeds  mb (earthquake) .  The  vertical  axis  is  the 
number  of  earthquakes/explosions  which  fall  into  a  given  magnitude 
difference  using  a  range  of  depths  (1-25  km),  source  orientations 
(all),  and  takeoff  angles  (10°-30°).  Surface  and  body  wave  mag¬ 
nitudes  were  averaged  over  azimuth. 
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We  find  the  following  relations  for  the  average  event  pair 
with  the  same  moment: 


and  in  general: 

1°9i0X  -  («;  -  Mj)-( mg  -  mj)  *  .35  .  (3.8) 

This  result  is  surprisingly  stable.  We  have  performed  this 
calculation  for  several  other  depth  ranges  (5-10  km,  15-20  km,  and 
25-50  km)  and  find  that  although  the  shape  of  the  distribution 
(Figure  5)  varies  somewhat,  the  peak  of  distribution  always  occurs 
in  a  narrow  range  witn  log^QX  *.35. 

For  small  events,  there  will  be  an  observational  bias  in  favor 
of  larger  amplitude  arrivals.  Small  amplitude  body  and  surface  wave 
arrivals  may  not  De  reported.  Both  body  and  surface  wave  magnitudes 
are  likely  to  oe  overestimated.  If  body  and  surface  waves  are 
observed  only  near  their  radiation  pattern  maxima,  then  the  source 
mechanism  effect  is  reduced  to  log^gXas.Z. 

Our  final  result  of  .35  for  the  magnitude  difference  is 
similar  to  Gilbert's  (1973)  estimate  of  .52,  but  the  interpretation 
is  quite  different.  He  suggested  that  for  a  given  moment,  body  wave 
magnitudes  should  be  almost  equal  while  surface  wave  magnitudes 
should  be  greater  for  earthquakes.  Gilbert's  analysis  depends  on 
approximating  the  double  dot  product  of  two  tensors  by  the  product 
of  their  Euclidian  norms.  Apparently,  this  procedure  does  not  yield 
good  separate  estimates  of  Rayleigh  and  P-wave  excitation  by 
earthquakes  relative  to  explosions  (i.e.,  A^/A*  and 

Afa/Ap)  ‘  th0u9h  it  does  give  a  reasonable  estimate  of  the 
relative  amplitude  factor  X. 
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FREQUENCY  DISTRIBUTION 


Figure 


.  Contour  plot  of  frequency  distribution  In  mb/Ms  plane.  Vertical 
axis  is  the  difference  between  mb  earthquakes  and  mb  of  explo¬ 
sions.  Horizontal  axis  Is  difference  between  Ms  of  earthquakes  and 
Ms  of  explosions.  In  general,  explosions  have  slightly  higher  Ms 
and  substantially  higher  mb  than  earthquakes  with  the  same  moment. 
The  frequency  of  occurrence  of  event  pairs  on  the  outermost 
contour  is  5  percent  of  the  maximum.  Inner  contours  are  multiples 
of  5  percent. 
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Figures  4  and  5  represent  our  best  estimate  of  the  difference 
between,  and  dispersion  of,  tne  earthquake  and  explosion  populations 
based  on  source  mechanism  alone.  Clearly,  source  mechanism  makes  a 
significant  contribution,  but  is  not  sufficient  to  explain  observed 
differences  between  m^  and  of  earthquakes  and  explosions. 
The  width  of  the  distribution  function  does  explain  why  there  should 
be  more  scatter  in  the  earthquake  population  and  why  earthquakes  may 
occasionally  look  like  explosions. 


IV.  EFFECT  OF  DEPTH 


The  depth  of  the  source  can  affect  both  body  and  surface  wave 
generation.  The  most  important  effect  of  depth  is  interference 
caused  Dy  reflected  phases.  The  pP  phase  from  an  explosion  is 

almost  identical  to  direct  P  but  opposite  in  sign  and  delayed  in 
time.  Tne  spectral  amplitude  is  cnanged  by  a  factor  of: 

j  1  +  A1  exp(-  (2i  u  h‘/a)  cose)  |  (4.1) 

where  A^  is  the  reflection  coefficient  of  the  free  surface  and 
Up  is  the  angular  frequency  of  the  P-wave.  For  steeply  descending 
rays  in  an  elastic  half-space  A^  2:  _l.  if  tne  spectrum  of  the 

direct  P-wave  is  g(u),  the  spectral  amplitude  of  P  +  pP  is  given 

approximately  by: 

|  9t(»)|  *  2] g(u») 1 1  sin  |  (4.2) 

where  t  is  tne  delay  time  xs=2h'/a.  At  low  frequencies  there  is 
complete  cancellation  while  at  frequencies  of  f  -  u/2u  -  (2n  +  l)/2x 
where  n  is  an  integer,  tne  amplitude  is  doubled  (Figure  6).  There 
will  be  cancellation  again  at  f  *  n/T.  Bache,  Day  and  Savino  (1979) 
found  time  delays  for  a  series  of  NTS  explosions,  with  mb  ranging 
from  5.5  to  6.5,  to  range  from  0.6  seconds  to  1.0  seconds.  A 

shallow  explosion  in  a  high  velocity  material  such  as  PILEDRIVER 
could  have  a  time  delay  as  short  as  0.2  seconds.  A  deep  explosion 
in  a  low  velocity  material  could  have  a  time  delay  as  long  as  2 
seconds.  Obviously  this  can  be  very  important  to  m^  measurements 
and  to  VFM  discrimination.  In  practice,  for  moderate  sized 
explosions,  interference  of  pP  is  more  likely  to  increase  m^  than 
to  decrease  it.  The  reason  is  that  the  amplitude  will  be  measured 
directly  from  the  seismogram  and  then  divided  by  the  instrument 
response  corresponding  to  the  measured  frequency.  This  frequency 
will  move  away  from  a  minimum  and  toward  a  local  maximum. 
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Figure  6.  Spectral  amplitude  of  direct  P  plus  elastic  pP  with  a  time 
delay  of  0.5  seconds.  The  two  waves  Interfere  destructively 
at  low  frequencies  and  alternately  interfere  constructively 
and  destructively  at  high  frequencies.  The  spectrum  of  the 
far  field  P-wave  is  the  spectrum  of  the  direct  P-wave  multi¬ 
plied  by  this  function. 
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The  effect  of  reflected  phases  on  earthquake  P-waves  is  more 
complex,  but  in  most  cases  less  Important.  If  the  earthquake  is 
more  than  a  few  kilometers  deep,  the  reflected  phases  will  be 
separated  in  time  from  the  direct  P-wave.  Because  of  the  radiation 
pattern  difference,  P  and  pP  will  usually  not  have  the  same 
amplitude.  Only  for  a  few  preferred  fault  orientations  will  there 
be  strong  destructive  or  constructive  interference. 

In  this  analysis  and  subsequent  examples,  we  use  a  linear, 
elastic  reflection  coefficient  to  determine  pP.  In  fact,  the 
upgoing  ana  reflected  P-waves  from  underground  explosions  travel 
through  a  highly  cracked  region  and,  in  most  cases,  cause  spallation 
of  the  upper  layers.  The  resulting  effect  on  far  field  body  waves 
is  not  wel 1 -determined  at  present.  Sache,  Day,  and  Savino  (1979)  in 
a  study  of  NTS  explosions  estimated  that  nonlinear  effects  reduce  pP 
by  a  factor  of  0.5  -  0.75.  On  the  other  hand,  Burdick  and 
Helmberger  (1979),  using  Russian  explosions,  found  that  the  body 
waves  could  be  modeled  using  a  linear  pP.  Day,  Rimer,  and  Cherry 
(1982)  show  that  conservation  of  momentum  requires  that  nonlinear 
effects  leave  pP  unchanged  at  long  periods.  The  most  likely  effect 
of  nonlinearities  is  attenuation  of  pP  at  high  frequencies.  This 
may  reduce  the  amount  of  interference  with  direct  P. 

For  small  explosions,  interference  of  pP  is  likely  to  reduce 

m.  .  Using  the  Mueller-Murphy  empirical  relation  for  depth 
D  1 

h  a  . 122W  xm,  the  travel  time  delay  is 

x  -  ,244W^3  seconds, 
a 

This  will  cause  destructive  interference  at  1  Hz  for  yields  of  50  KT 
or  less  in  granite  and  for  5  KT  or  less  in  tuff. 

The  effect  of  pP  interference  is  to  increase  the  separation  of 
the  earthquake  and  explosion  populations  on  the  plane  for 

moderate  and  high  yields  by  as  much  as  .3,  and  to  reduce  the 


27 


S-CUBCD 


1 


separation  at  low  yields.  This  causes  an  apparent  convergence  of 
the  populations  at  small  magnitudes  which  will  be  evident  in  the 
synthetic  seismogram  results  presented  in  a  later  section. 

The  effect  of  reflected  phases  on  the  variable  frequency 
magnitude  discriminant  is  the  opposite  of  the  effect  on  the 
mb:Ms  discriminant.  A  summary  of  the  VFM  method  is  given  by 
Bache,  Day  and  Savino  (1979).  The  method  consists  of  applying  a 
"comb"  of  narrow  band  filters  over  a  range  of  frequencies  from 
0.4  Hz  to  4.0  Hz  to  a  snort-period  seismogram.  Magnitude 
measurements  are  then  made  at  two  frequencies  as  far  apart  as  noise 
will  permit  —  usually  about  0.5  Hz  and  3.0  Hz.  In  practice,  the 
spectra  are  smoothed  to  remove  some  of  the  dips  in  the  spectra. 
This  nas  tne  result  that  interference  of  pP  tends  to  amplify  the 
spectra,  especially  at  higher  frequencies.  If  the  first 

interference  peak  occurs  at  0.5  Hz  or  less,  the  effects  should 
average  out  and  cause  no  apparent  amplification.  On  the  other  hand, 
for  small  events,  destructive  interference  at  0.5  Hz  will  reduce  the 
low  frequency  amplitude  for  shallow  explosions  with  pP  -  P  time 
delays  less  than  about  0.3  seconds,  increasing  the  separation  by  as 
much  as  0.5.  Thus,  the  interference  of  pP  will  aid  in  the  explosion 
identification  using  the  VFM  discriminant,  and  hinder  the 
identification  using  mb:Ms.  For  deeper  events,  however, 
interference  of  pP  will  help  in  the  identification  of  explosions  by 
mb:Ms  and  complicate  the  identification  (by  adding  spectral  dips 
ana  peaks)  using  the  VFM  discriminant. 

Surface  wave  generation  is  also  affected  by  source  depth.  A 
vertical  dip-slip  earthquake  at  the  surface  causes  no  surface  wave 
generation.  A  strike  slip  or  45  degree  thrust  earthquake  at  a 
critical  depth,  or  observed  at  a  critical  frequency  given 
approximately  oy  wh  as  cR,  also  causes  no  surface  wave  generation 
(This  is  true  for  a  layered  half-space  as  well  as  for  a  uniform 
half-space  although  the  critical  depths  will  be  different.  See 
Tade  2).  In  Figure  7  we  have  plotted  the  relative  magnitude  factor 
(from  Equation  3.5)  as  a  function  of  frequency  for  strike  slip  and 
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LOG  ( A^/Aj) 

Depth  =  10  km  Depth  =  I  km 


FREQUENCY  DEPENDENCE  OF  RAYLEIGH  GENERATION 


Rayleigh  Wave  Frequency 


Figure  7.  Relative  Rayleigh  wave  amplitudes  as  a  function  of  frequency,  depth 
and  dip  angle.  Observation  point  for  dip-slip  is  90°,  for  strike- 
slip  45°  close  to  maxima  in  the  radiation  pattern.  Takeoff  angle 
for  body  waves  is  10°  (body  waves  are  not  frequency  dependent).  The 
Rayleigh  wave  velocity  Is  3.185  km/sec.  Vertical  axis  is  1og(A^/A£) 
for  same  body  wave  amplitude. 
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dip-slip  earthquakes  at  two  depths,  for  a  number  of  dip  angles. 
There  is  a  clear  frequency  and  depth  dependence  to  the  Rayleigh  wave 
amplitude.  Tsai  and  Aki  (1971)  suggested  that  this  depth/ frequency 
dependence  could  be  used  to  discriminate  earthquakes  from 
explosions.  Unfortunately,  as  Figure  7  shows,  the  frequency 
dependence  is  different  for  different  depths  and  different 
double-couple  orientations.  In  particular,  the  dips  in  the  spectrum 
are  only  obvious  for  source  orientations  close  to  vertical  strike 
slip,  vertical  dip-slip,  or  45  degree  dip-slip.  The  spectrum  is 
relatively  smooth  and  homogeneous  at  intermediate  orientations. 

It  is  customary  in  source  studies  to  consider  these  three 
source  orientations  individually  and  then  to  make  conclusions  based 
on  the  results  using  the  fact  that  all  double-couple  sources  are 
linear  combinations  of  these.  As  we  have  seen,  however,  these 
particular  source  orientations  have  some  anomalous  features  which, 
in  general,  do  not  occur  for  linear  combinations  of  the  sources;  so 
care  is  needed  in  making  general  conclusions  from  properties  of 
these  sources  alone. 
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V.  EFFECT  OF  SOURCE  REGION  ELASTIC  PROPERTIES 


We  have  shown  that  source  mechanism  and  source  depth  can 
explain  part,  but  not  all,  of  the  observed  mb-Ms  difference 
between  earthquakes  and  explosions.  Source  depth  can  also  explain 
part  of  the  difference  fount,  using  the  VFM  discriminant.  The 

remaining  differences  may  be  explained  by  differences  in  the  source 
spectra  and  the  material  properties  of  the  source  region.  In  this 
section  we  consider  the  effects  of  the  latter. 

Actually  the  near-surface  geologic  structure  has  two  very 
different  influences  on  the  ixplosion  signal  which  we  must 
distinguish.  First,  the  nonlinear  rheology  of  the  medium  affects 

the  explosion  source  spectrum.  This  effect  is  considered  in  the 
next  section.  In  this  section,  we  are  concerned  with  the  manner  In 
wnich  near-source  elastic  properties  modify  transmission  of  the 
explosion  signal.  Here,  we  take  a  simplified  view  of  the  wave 

propagation  in  order  to  estimate  this  transmission  effect  alone.  In 
a  later  section,  we  compute  synthetic  seismograms  which  incorporate 
the  effect  of  crustal  geology  rigorously. 

We  model  the  source  region  as  a  layer  of  lower  velocity 

material  over  a  half-space.  We  restrict  attention  to  explosions  and 
earthquakes  occurring  in  the  same  general  area;  so  we  will  not 

concern  ourselves  with  the  travel  path  outside  of  the  source 
region.  Since  eartnquakes  are,  in  general,  deeper  than  explosions, 
we  will  consider  the  earthquakes  to  be  located  in  the  half- space 
while  explosions  will  be  located  in  the  overlying  layer. 

We  need  to  estimate  the  body  wave  and  surface  wave  amplitudes 
for  an  explosion  or  earthquake  in  this  two  layered  structure.  We 

will  consider  only  the  direct  (far-field)  P-wave  since  reflected 
phases  have  been  considered  separately.  For  the  earthquake,  the 
body  wave  will  be  equal  to  the  free  space  body  wave  for  which  the 
spectral  amplitude  a£  Is  given  by  Equation  3.1.  For  the 

explosion,  however,  the  free-field  body  wave.  Equation  3.2,  is 


modified  by  tne  transmission  coefficient  from  the  upper  to  the  lower 
layer.  For  steeply  descending  rays,  the  transmission  coefficient 
(T  )  is  approximately 

,  2 

T  2ttx  °x 

xq  *  ®q  <Vx"  *  ®qV 

where  0  and  a,  are  the  density  and  P-wave  velocity  of  the 

q  q 

half-space  and  and  ax  are  the  density  and  P-wave  velocity  of 
the  layer  (Cerveny  and  Ravindra,  1971).  The  continuous  limit 
(gradual  transition  between  materials)  may  also  be  found  using  WKBJ 
theory  (Aki  and  Richards,  Chapter  9,  page  *24).  The  body  wave 
spectrum  from  a  source  with  spectrum  mx(w)  is  given  for  a  steeply 
descending  ray  by 


m  (M) 

u(«)  -  - 5 - m -  (5.1) 

4*  taq  ®q  °x  A)  * 

where  R  is  a  geometrical  spreading  factor,  [n  our  case. 


The  amplitude  j  u(w)  |  is  the  explosion  body  wave  amplitude  a£. 
Comparing  with  Equation  3.2,  we  find  that  the  transmission 
coefficient  is 


These  two  transmission  coefficients  are  similar  In  amplitude, 
especially  for  velocity  differences  of  less  than  a  factor  of  two. 
4e  will  use  the  continuous  coefficient  because  of  Its  simpler  form. 
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Aki  and  Ricnards  (1980,  Chapter  7)  give  a  general  expression 
for  the  excitation  of  surface  waves  by  a  source  with  a  given  seismic 
moment  tensor  In  terms  of  the  free  eigenfunctions  of  a  layered 
medium.  For  the  case  of  an  explosion,  the  Rayleigh  wave  amplitude 
is  given  0y 


ux  -  F(r) 


Kr^h) 


r4(h) 

*  — r~ 


mx(w) 


(5.3) 


where  F(r)  is  a  propagation  function,  K  Is  the  Rayleigh  wave  number, 
a,  s,  and  p  are  the  P-wave  and  $-wave  velocities  and  density  of  the 
source  material,  h  is  tne  source  depth  and  r^  and  r4  are, 
respectively,  horizontal  displacement  and  normal  stress 
eigenfunctions  of  the  layered  space.  The  stress  eigenfunction  r4 
vanishes  at  tne  free  surface.  For  periods  at  which  M$  is  measured 
(about  20  seconds),  r4  is  small  for  all  normal  explosion  depths. 
Furthermore,  at  such  periods,  r^  is  nearly  equal  to  the  half-space 
eigenfunction.  Also,  we  can  neglect  tne  effect  of  source  material 
on  explosion  source  spectral  shape  at  these  long  periods.  The 
amplitude  of  surface  waves  generated  by  explosions  in  different 
materials  will  therefore  be  proportional  to  /<*  .  For 

different  media  with  similar  Poisson's  ratios,  the  amplitude  of 
surface  waves  generated  by  explosions  will  be  proportional  to  the 
explosion  moment  M*. 

Combining  the  results  for  body  and  surface  waves,  we  find  that 
the  contribution  of  source-region  material  properties  alone  to  the 
relative  amplitude  factor  is  given  by 


(5.4) 


The  logarithm  of  this  number  is  the  contribution  of  source-region 
elastic  properties  alone  to  tne  distance  between  the  earthquake  and 
explosion  populations  on  an  m^:.^  diagram.  It  is  added  to  the 
separation  found  before  due  to  radiation  pattern  which  was  estimated 


33 


SCUBED 


I 


earlier.  This  contribution  can  be  very  significant.  For  explosions 
in  low  velocity  materials  such  as  tuff,  it  can  increase  the 
magnitude  difference  by  about  .6,  assuming  that  earthquakes 
generally  occur  in  the  crystalline  basement. 
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VI.  SPECTRAL  COMPARISON  OF  SOURCES 

f  In  this  section  we  will  make  synthetic  m^Mg  curves  and 

synthetic  VFM  curves  Oy  combining  the  contributions  analyzed  in  the 
previous  sections.  That  analysis  is  sufficient  to  give  estimates  of 
relative  spectral  amplitude.  The  first  step  in  the  synthesis  is  to 
*  turn  tnese  relative  spectral  amplitudes  into  estimates  of  body  and 
surface  wave  magnitudes.  To  do  this  we  need  at  least  one  reference 
point.  We  use  the  result  of  Murphy  (1977),  based  on  NTS  data,  that 
a  100  tCT  explosion  in  a  tuff /rhyolite  medium  has  an  of  5.5  and 
an  M$  of  4.0.  Using  these  values  and  the  spectrum  for  the 

Mueller-Murpny  Tuff/Rhyolite  explosion  source,  we  can  scale  all  of 
the  explosions: 

*2  -  ]°ho  +  c£  <a> 

(6.1) 

M;  -  1og10  (A*)  ♦  C*  (b) 

where  the  constants  c£  and  C*  are  chosen  to  give  agreement 
with  Murphy's  result  for  100  KT.  In  this  section,  we  will 
approximate  a£  using  Equation  5.1  so  that  both  the  source 
spectrum  and  the  source  region  transmission  properties  are  accounted 
for.  As  discussed  In  the  last  section,  we  can  simply  take  A* 
proportional  to  the  source  spectrum  mx< 

For  earthquakes,  we  use  expressions  analogous  to  Equation  6.1, 
out  modify  the  constants  to  include  the  effects  of  source  mechanism 
and  radiation  pattern.  For  mjj  and  M*.  we  write 


mb  ■ 

lo,10  AS*cb 

(a) 

mJ  . 

>°S10  a2  *  $ 

(b) 
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and  from  our  earlier  analysis  we  know  that  c£  **  eg  +  .45 

and  C*  «  +  .1.  We  will  evaluate  Ag  and  Ag 

t  from  the  theoretical  RMS  earthquake  spectra. 

The  functions  Afa  are  evaluated  at  a  period  of  one  second 
while  As  are  evaluated  at  20  seconds.  These  amplitudes  will,  of 
course,  vary  from  the  measured  time-domain  amplitude  because  of  the 
*  finite  spectral  width  of  the  seismometer  response,  but  we  will 

assume  for  now  that  this  is  an  adequate  approximation.  In  the  next 
section  we  will  use  synthetic  seismograms  to  test  this  approximation. 

Using  MKS  units  for  all  quantities  and  using  unit  distance  for 
the  boay  waves,  we  find  the  constants  in  Equations  6.1  and  6.2  to  be 

c£  *  4.75 

eg  *  4.3 

eg  *-12.1 

Cg  *-12.2  ; 

since,  in  the  long  period  limit,  the  surface  wave  amplitudes  are 
proportional  to  the  seismic  moment,  this  also  gives  a 
moment/magnitude  relation.  For  small  explosions 

Ms  *  log10  “  12a  (6-3) 

and  for  small  earthquakes: 

«J  *logiQ  -  12.2  .  (6.4) 

We  can  compare  MS:MQ  relations  obtained  from  Equations  6.1 
and  6.2  with  empirical  observations  of  Ms  and  moment.  Figure  8 
shows  measurements  for  a  large  number  of  earthquakes  using  data  from 
Kanamori  and  Anderson  (1975)  and  Hanks  and  Kanamori  (1979).  The 
solid  curves  on  the  figures  are  the  theoretical  values  of 
computed  from  Equation  6.2  using  the  6  finite  difference  earthquake 
1  calculations.  In  general,  agreement  with  the  earthquake 

observations  is  quite  good  in  spite  of  the  fact  that  we  normalized 
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MS 


Figure  8.  Moment  versus  surface  wave  magnitude.  Ms  is  computed  from  the 
scaled  20  second  amplitude  for  the  six  earthquake  spectra.  Oata 
points  are  from  Kanamori  and  Anderson,  1975  (-*-)  and  Hanks  and 
Kanamori,  1979  (x). 
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our  magnitudes  to  an  explosion  model.  Geller  (1976),  using  a 
Haskell  dislocation  model  and  the  same  data  set,  estimated 
=  -11.9. 

A  simple  Haskell  dislocation  model  predicts  a  slope  of  one  at 
small  magnitudes  followed  by  a  slope  of  3/2  at  intermediate 
magnitudes,  and,  eventually,  an  infinite  slope  on  the 
moment/magnitude  plot.  Kanamori  and  Anderson  (1975)  determined  that 
tne  slope  should  be  3/2  for  most  events  with  M  between  6  and  8, 
and  used  this  relation  to  estimate  stress  drops.  Our  results  do  not 
show  a  clear  slope  of  3/2  at  any  point.  In  fact,  the  slopes  are 
nearly  one  for  Ms  £  7  for  all  theoretical  curves  and  become  very 
steep  in  an  irregular  manner  for  larger  earthquakes.  For  <  7, 
tnerefore,  the  curves  are  very  insensitive  to  changes  in  stress  drop. 

We  can  now  compute  synthetic  mt):Ms  curves  using  Equations 
6.1  and  6.2  for  all  of  the  explosion  and  earthquake  models.  The 
result  is  shown  in  Figure  9.  These  synthetic  curves  include  effects 
of  source  spectrum,  source  region  transmission  properties,  and 
average  source  mechanism,  but  include  no  correction  for  interference 
by  surface-reflected  phases.  The  solid  lines  on  the  figure  are  the 
earthquake  models;  the  dashed  lines  are  the  finite  difference 
explosion  models;  and  the  dotted  lines  are  Mueller-Murphy  explosion 
models.  The  explosion  population  is  very  well  separated  from  the 
earthquake  population  even  at  small  magnitudes.  There  is,  of 
course,  much  more  scatter  in  the  earthquake  population  than  is  shown 
since  we  have  used  only  the  peak  of  the  source  mechanism  frequency 
distribution  (Figure  4).  An  interesting  observation  is  that  for 
most  sources,  spectral  differences  are  of  secondary  importance  for 
events  with  <  4. 

As  mentioned  earlier,  fault  mechanism  is  responsible  for  an 
average  difference  of  about  .35  between  the  populations.  At  smail 
magnitudes,  most  of  the  remaining  difference  is  due  to  variations  in 
source  material  properties  as  discussed  earlier.  The  surface  to 
body  wave  amplitude  ratio  A$/Ab  is  proportional  to  so 

an  explosion  in  a  lower  velocity  material  has  M$  decreased 
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MB/MS  COMPARISON 


Figure  9 


MS 


comparison  using  source  spectra  at  one  and  twenty  seconds. 
Earthquakes  (solid  lines)  are  identified  by  number  as  in  Figure  1. 
Muel 1 er-Murphy  explosion  sources  (dotted  lines)  are  for  Salt  (S), 
Granite  (G),  Tuff-Rhyolite  (TR)  and  Shale  (H).  Numerical  explosion 
sources  (dashed  lines)  are  for  Salt  (S),  Dry  Granite  (GD),  Wet 
Granite  (GW)  and  Tuff  (T).  The  theoretical  curves  compare  well 
with  the  data  of  Marshall  and  Basham  (1972).  Separation  between 
the  populations  at  low  magnitude  is  due  primarily  to  source 
mechanism  and  source  region  elastic  properties.  This  figure 
does  not  include  interference  of  pP  and  uses  only  the  peak 
source  mechanism  effect.  The  distribution  in  Figure  4  must  be 
superimposed  on  the  earthquake  curves. 
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relative  to  m^.  This  accounts  for  an  average  magnitude  difference 
of  0.4  -  1.0  even  in  the  absence  of  spectral  differences.  Figure  8 
is  in  very  good  agreement  witn  the  empirical  data  of  Marshall  and 
Basham  (1972,  Figures  6  and  8). 

Although  source  spectral  differences  may  be  of  secondary 

importance  to  mtj:Ms  discrimination,  they  are  of  primary 
importance  to  VFM  discrimination.  Figures  10  to  12  are  comparisons 
of  spectra  of  individual  earthquake  models  and  explosion  models. 

The  explosion  models  were  all  evaluated  at  100  KT,  and  earthquake 
models  were  scaled  to  have  the  same  P-wave  spectral  amplitude,  in 

the  long  period  limit,  as  the  100  KT  explosion.  For  the  purpose  of 
these  comparisons,  radiation  pattern  effects  were  neglected 
(c£*c£),  and  the  maximum  earthquake  P-wave  amplitude  is 
used.  The  earthquake  spectral  amplitudes  are  root  mean  square 

spectral  amplitudes  normalized  to  the  maximum  of  the  radiation 
pattern.  The  solid  line  on  each  figure  is  the  logarithm  of  the 
P-wave  spectral  ratio  (explosion  to  earthquake).  The  normalization 
requires  that  this  ratio  equal  one  in  the  long  period  limit.  Also 
shown  on  the  figures  are  the  individual  explosion  (dashed)  and 
earthquake  (dotted)  spectra  normalized  to  one  at  long  period. 

These  figures  show  clearly  why  VFM  works  as  a  discriminant. 
In  all  cases  tne  explosions  have  significantly  more  high  frequency 
energy  than  tne  earthquakes.  In  the  first  set  of  figures.  Figure  10 
shows  the  Mueller-Murphy  tuff/rhyolite  explosion  spectrum  (at  100 
KT)  compared  to  each  of  the  earthquake  sources  while  Figure  11  shows 
earthquake  source  number  1,  the  spontaneous  rupture  with  a  100  bar 
stress-drop  in  a  circular  region  surrounded  by  a  zero  stress-drop 
region,  compared  with  the  four  numerical  explosion  sources.  Figure 
12  shows  earthquake  number  1  compared  with  the  four  empirical 
(Mueller-Murphy)  explosion  sources. 

The  earthquake  and  numerical  explosion  sources  may  be  scaled 
directly  to  otner  yields  by  simply  translating  the  axis  by 
1/3  logio(W/ 100) .  Scaling  for  the  Mueller-Murphy  sources  is  more 
complicated  (see  Appendix  A),  and  these  sources  become  more 
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Figure  10.  Spectral  comparison  of  earthquake  sources  with  alOO  KT  explosion 
in  Tuff  (Mueller-Murphy).  Dotted  line  is  earthquake  RMS  spec¬ 
trum  observed  at  a  radiation  pattern  maximum  and  normalized  to 
same  long  period  body  wave  limit  as  the  explosion  body  wave. 
Dashed  line  is  the  explosion  P-wave  spectrum;  solid  line  is  the 
logarithmic  difference  between  the  spectra.  In  all  cases,  the 
explosion  has  more  high-frequency  energy  than  the  earthquakes. 
Top  three  figures  correspond  to  the  first  three  earthquakes  in 
Figure  1,  and  bottom  three  to  the  last  three. 
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Figure  11.  Spectral  comparison  of  numerical  explosion  models  at  100  KT  with 
earthquake  number  1  (approximately  circular  with  100  bar  stress 
concentration).  Again,  all  explosions  have  more  high-frequency 
»  energy  than  the  earthquake  although  some  complexity  is  observed 

at  very  high-frequency. 
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Figure  12.  Spectral  comparison  of  Muell er-Murphy  explosion  models  at 
100  KT  with  earthquake  number  1. 
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earthquake-like  at  smaller  yields.  Figure  13  shows  the 

Mueller-Murphy  explosion  source  for  tuff/rhyolite  at  a  yield  of 
2.5  KT  compared  with  earthquake  number  3.  On  the  basis  of  source 
spectrum,  it  would  appear  that  very  small  explosions  in  low  velocity 
materials  may  not  oe  discriminated  by  VFM.  However,  at  the  shallow 
depth  (165  meters)  assumed  for  this  source,  interference  with  pP 
would  allow  discrimination. 

The  next  set  of  figures,  Figures  14  to  16,  are  synthetic 
variable  frequency  magnitude  curves  made  by  evaluating  Equations 
6.1a  and  6.2a  at  two  frequencies  and  scaling  to  a  number  of 

different  amplitudes.  Cube  root  scaling,  as  described  earlier,  was 

used  for  the  numerical  earthquake  and  explosion  sources.  The 

Mueller-Murphy  sources  were  scaled  with  yield  as  described  in 

Appendix  A.  The  effect  of  radiation  pattern  was  included  by 
reducing  mb  for  earthquakes  by  using  the  c£  and  c£  given 
before  Equation  6.3, 

The  first  pair  of  figures  (Figure  14)  are  for  a  low  frequency 

level  f^  »  .05  Hz  and  high  frequency  levels  of  f^  «  1  Hz  and 

2  Hz.  The  earthquake  and  explosion  populations  are  separated  on 

both  plots,  but  the  separation  is  quite  small  for  <  4.5  with 
f2  -  1  Hz.  Using  f2  »  2  Hz,  the  separation  is  much  better.  The 

populations  are  clearly  separated  for  >  4.0. 

In  practice,  because  the  noise  level  increases  at  low 
frequencies,  it  is  not  possible  to  use  frequencies  much  less  than  .5 
Hz.  Figure  15  shows  synthetic  VFM  curves  for  f^  -  .5  Hz  and 

f2  »  2,  3  and  4  Hz.  For  small  events,  separation  is  almost  as 

good  as  for  f^  *  .05  Hz.  For  f2  ■  4  Hz,  separation  is  quite 

clear  even  down  to  m^  »  3.5.  For  very  large  events,  separation  is 
not  as  good.  This  is  because  the  low  frequency  level  f^  is  moving 
above  the  corner  frequency  and  seeing  the  high  frequency  complexity 
of  the  spectrum.  Separation  is  especially  good  for  explosions  in 
higher  velocity  materials.  Explosions  in  granite  are  much  more 
clearly  separated  than  explosions  in  tuff. 


44 


S-COBED 


dWU  91)  I 


_ 

I 


I 


I 


LOG  FR£G 


* 

Figure  13.  Spectral  comparison  of  2.5  KT  Mueller-Murphy  Tuff  explosion  with 
earthquake  number  3.  The  only  case  In  which  explosions  do  not 
clearly  have  more  high  frequency  energy  than  earthquakes  Is  for 
the  empirical  model  of  small  explosions  in  low  velocity  materials. 
Again,  the  dotted  line  is  the  earthquake,  the  dashed  line  is  the 
f  explosion,  and  the  solid  line  Is  the  difference. 
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Figure  14.  Synthetic  variable  frequency  magnitude  curves  using  f.  =0.05  Hz  for  the  low  frequency  and 
f2  *  1.0  Hz  and  2.0  Hz  for  the  high  frequency.  Sol  id' lines  are  earthquakes  marked  at 
moment  of  10'7  NT-M.  Dashed  lines  are  numerical  explosion  sources;  dotted  are  Mueller- 
Murphy  explosion  sources  marked  at  100  KT  yield.  The  anomalously  low  amplitude  mark  is 
for  dry  granite.  Separation  of  the  earthquake  and  explosion  populations  is  much  better 
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Figure  15.  Synthetic  Variable  Frequency  Magnitude  curves  for  f.  -  0.5  Hz,  f?  =  2.0,  3.0,  and  4.0  Hz. 

These  are  commonly  used  frequencies.  The  noise  levil  increases  rapidly  below  0.5  Hz,  and 
observations  become  difficult  above  4.0  Hz. 


frequency  moves  below  0.8  Hz . 


It  could  oe  advantageous  to  use  a  higher  low-frequency  level 
to  avoid  low-frequency  noise.  Figure  16  shows  synthetic  VFM  curves 
for  f^  -  .8  Hz  and  ■  2,  3  and  4  Hz.  The  earthquake  and 

explosion  populations  are  still  fairly  well  separated,  but  the 

region  of  clear  separation  is  smaller  than  was  the  case  with 
f^  -  .5  Hz.  VFM  measurements  using  f^  a  .8  Hz  would  be  able  to 
separate  most  events,  but  there  is  a  greater  chance  of 
misidentification  tnan  there  is  using  lower  values  of  f^. 

The  synthetic  VFM  curves  are  quite  insensitive  to  the 

assumptions  we  have  made.  A  change  in  stress-drop  of  the 

earthquakes,  for  example,  simply  moves  the  earthquake  curves  along  a 
45  degree  angle  on  the  plots.  Stress  drops  of  1  -  5  Kb  are 
generally  required  to  cause  mixing  of  the  earthquake  and  explosion 
populations.  Explosions  in  granite  only  mix  with  very  large,  very 
high  stress '  drop  earthquakes.  The  numerical  granite  explosion 
models  only  mix  with  earthquakes  with  stress  drops  greater  than 
10  Kb  and  mb  >  5.5.  Explosions  in  tuff  do  not  discriminate  as 
easily.  There  is  some  mixing  with  a  stress  drop  of  250  bars  at 
mb  «  4.5.  Similarly,  radiation  pattern  effects  shift  the 

earthquakes  up  or  down  along  a  45  degree  angle.  The  curves  shown 
are  for  the  average  of  the  radiation  pattern.  At  the  maximum,  there 
is  a  small  amount  of  mixing  with  the  Muel ler-Murphy  tuff  source  at 
low  yields  crossing  into  part  of  the  earthquake  population  as  the 
spectral  comparison  of  Figure  13  showed.  Earthquake  Models  1,  3, 
and  4  m  x  with  the  explosions  at  much  lower  stress  drop  than 
Earthquake  Models  2,  5,  and  6. 

We  have  used  spectra  evaluated  exactly  at  the  observation 
frequency  for  this  comparison.  In  practice,  any  measurement  is 
taken  over  a  finite  frequency  window  and  will  not  be  centered 
exactly  at  the  desired  frequency.  VFM  studies,  however,  use  a  very 
narrow  band  width;  so  this  should  be  quite  a  good  approximation. 

Probably  the  most  important  approximation  made  here  was  the 
neglect  of  reflected  body  wave  phases.  As  discussed  earlier,  these 
are  expected  to  preferential ly  suppress  ( ^^ ) »  the  low 

frequency  magnitude,  which  will  enhance  event  separation. 


VII.  TIME -DOMAIN  SOURCE  COMPARISON 


In  this  section  we  make  synthetic  seismograms  from  our 

earthquake  and  explosion  source  functions  and  make  measurements  of 
mb  in  the  time-domain.  While  most  of  the  essential  source 

characteristics  are  contained  in  the  preceding  spectral  comparisons, 
spectral  estimates  of  m^  will  not,  in  general,  be  equal  to 

time-domain  measurements  of  m^.  Time-domain  amplitudes  are 
measured  using  an  instrument  with  a  f inite-frequency  band  width  so 
the  measured  frequency  will  be  shifted  away  from  the  dominant 
instrument  frequency  in  the  direction  of  increasing  spectral 
amplitude. 

Synthetic  seismograms  include  the  effects  of  source  radiation 
pattern  and  interference  of  reflected  and  transmitted  phases  which 
were  neglected  or  approximated  earlier.  The  spectral  comparisons  in 
the  last  section  used  the  root  mean  square  spectra  for  the 
earthquakes,  which  lack  some  of  the  complexity  of  the 

position-dependent  individual  spectra.  Time-domain  comparisons 
include  tne  phase  information  which  was  not  used  in  the  spectral 
comparisons. 

To  see  how  time-domain  measurements  compare  with 
frequency-domain  measurements,  we  made  synthetic  seismograms  for  a 
rwmoer  of  the  earthquake  and  explosion  sources.  The  seismograms 
were  constructed  using  the  modification  of  the  method  of  Bache  and 
Harkrider  (1976)  described  in  Appendix  B. 

The  following  parameters  were  used  for  all  the  earthquake 
syntnetics.  The  earthquake  is  at  a  depth  of  5  kilometers  in  a 
structure  similar  to  the  eastern  U.S.  structure  SI  determined  by 

Bache,  Swanger  and  Shkoller,  1980,  which  is  listed  in  Table  3.  The 

receiver  was  located  in  a  4  km/sec  half-space  at  a  distance  of  35 
degrees  (3900  km);  so  the  takeoff  angle  from  the  source  was  27 
degrees.  A  t*  of  .8  was  used  for  attenuation.  The  fault  was  placed 
in  45  degree  dip-slip  orientation  (strike  -  0*,  dip  »  45*,  rake  » 
90*)  and  observed  at  an  azimuth  of  90  degrees,  near  a  maximum  in  the 
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SOURCE  REGION  STRUCTURE  FOR  SYNTHETICS 


I 

Layer 

Depth  (km) 

a 

(km/sec) 

6 

(km/ sec) 

(q/cnr) 

1 

0.6 

3.7 

2.2 

2.1 

2 

2.6 

4.5 

2.5 

2.2 

i 

3 

4.1 

5.6 

3.2 

2.6 

4 

6.2 

6.0 

3.46 

2.7 

5 

13.2 

6.3 

3.6 

2.8 

6 

20.0 

6.4 

3.7 

2.9 

7 

oc 

6.6 

3.8 

3.0 

I 


51 


I 


SOUSED 


radiation  pattern.  A  few  synthetics  were  also  made  at  an  azimuth  of 
270  degrees.  Since  we  are  interested  primarily  in  differences  in 
the  source  and  source  region,  a  geometric  spreading  factor  of  1.2  x 
10  was  used  in  lieu  of  a  mantle  transfer  function.  A  WWSSN 
short  period  instrument  response  was  used.  The  maximum  response  of 
this  instrument  occurs  at  a  period  of  .62  seconds. 

We  also  made  synthetic  seismograms  for  the  Mueller-Murphy 
explosion  source  functions  in  four  different  materials.  The 
synthetic  seismograms  were  made  in  the  same  manner  as  the  earthquake 
seismograms  using  the  method  of  Sache  and  Harkrider,  except  that  the 
source  region  material  constants  were  modified.  The  top  two  layers, 
to  a  depth  of  2.6  km,  were  taken  to  have  the  same  material 

properties  as  the  explosion  source  material  (as  in  Table  1).  The 

1/3 

Mueller-Murphy  depth  of  h  a  122  W  ,  where  h  is  in  meters  and  W 
is  in  kilotons,  was  used  for  the  explosion  source  depth. 

Magnitude  measurements  were  made  to  resemble,  as  much  as 
possiDle,  magnitude  measurements  from  observed  seismograms.  For 
most  cases,  the  amplitude  used  was  the  peak  to  peak  amplitude  from 
the  first  extremum  to  the  second  extremum.  The  period  used  was 
twice  tne  time  interval  between  these  points.  An  exception  to  this 
metnod  of  measurement  was  made  for  some  of  the  larger  earthquakes. 
For  these  events  the  pulse  duration  was  long  enough  that  the 
amplitude  became  much  larger  after  the  first  cycle.  In  these  cases, 
the  largest  amplitude  among  the  first  few  cycles  was  used.  At  the 
distance  (35  degrees)  used  for  the  synthetics,  the  body  wave 
magnitude  mb  is  defined  by 

mb  *  lo9lo(£r)  +  3*4 

where  G  is  the  gain  of  the  instrument  at  period  T,  and  A  is  the  peak 
to  peak  amplitude  in  nanometers. 

Measurements  made  from  tne  synthetic  seismograms  are 
summarized  in  Tables  4  (earthquakes)  and  5  (explosions).  Shown  on 
the  tables  are  tne  moment  of  the  event,  m^  estimated  from  the 
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TABLE  4 


MAGNITUDE  ESTIMATES  FROM  SYNTHETIC  EARTHQUAKE  SEISMOGRAMS 


Event  # 

1og10(Mo) 

mb(f) 

T/n 

"b<T> 

1 

17.10 

5.11 

1.04/1 

5.76 

17.10 

5.11 

1.12/lb 

5.73 

17.00 

5.06 

1.00/1 

5.68 

16.40 

4.66 

.91/1 

5.21 

15.80 

4.15 

.86/1 

4.68 

14.50 

2.87 

.82/1 

3.43 

2 

18.20 

5.13 

1.30/1 

6.24 

17.38 

5.00 

1.13/1 

5.78 

17.38 

5.00 

1.17/lb 

5.77 

16.60 

4.64 

1.02/1 

5.27 

16.00 

4.24 

.92/1 

4.80 

15.00 

3.38 

.84/1 

3.91 

3 

18.45 

5.67 

1.87/1 

6.80 

18.45 

5.67 

1.62/lb 

6.64 

17.70 

5.44 

1.24/1 

6.32 

16.90 

5.09 

.94/1 

5.70 

16.30 

4.62 

.87/1 

5.18 

15.70 

4.08 

.83/1 

4.62 

4. 

18.00 

5.54 

1.64/1 

6.42 

18.00 

5.54 

1.12/lb 

6.32 

17.50 

5.30 

1.38/1 

6.12 

16.80 

4.95 

.99/1 

5.56 

16.20 

4.50 

.89/1 

5.06 

15.60 

3.97 

.84/1 

4.52 

53 


S CUBED 


TABLE  4  (continued) 

MAGNITUDE  ESTIMATES  FROM  SYNTHETIC  EARTHQUAKE  SEISMOGRAMS 


Event  # 


5 


log10(Mo) 

mb(f) 

T/n 

VT> 

18.10 

5.23 

1.28/5 

6.52 

17.15 

4.90 

1.07/1 

5.38 

17.15 

4.90 

1.16/lb 

5.64 

16.70 

4.68 

1.21/1 

5.20 

16.10 

4.25 

1.04/1 

4.82 

15.40 

3.70 

.89/1 

4.25 

18.35 

5.39 

1.68/4 

6.55 

18.35 

5.39 

1-1 3/ 4b 

6.21 

17.80 

4.84 

1.78/1 

6.09 

16.50 

4.59 

1.02/1 

5.22 

15.90 

4.16 

.89/1 

4.74 

15.30 

3.65 

.85/1 

4.20 
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TABLE  5 

MAGNITUDE  ESTIMATES  FROM  SYNTHETIC  EXPLOSION  SEISMOGRAMS 
Material  Yield  mb^  T 


Salt 

1 

14.91 

4.28 

.71 

4.07 

10 

15.68 

5.15 

.73 

5.09 

100 

16.44 

5.94 

.77 

6.07 

1000 

17.21 

6.58 

.97 

6.93 

Tuff/ Rhyolite  1 

14.59 

3.92 

.70 

3.76 

10 

15.35 

4.70 

.72 

4.79 

100 

16.10 

5.50 

.82 

5.74 

1000 

16.87 

6.24 

1.06 

6.53 

Granite 

1 

14.87 

3.88 

.75 

3.72 

10 

15.63 

4.72 

.71 

4.73 

100 

16.39 

5.65 

.76 

5.78 

1000 

17.16 

6.48 

.89 

6.75 

Shale 

1 

14.64 

3.81 

.72 

3.63 

10 

15.40 

4.62 

.72 

4.67 

100 

16.16 

5.49 

.76 

5.69 

1000 

16.92 

6.28 

.97 

6.59 
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frequency-domain  measurements  in  the  last  section  (at  one  Hertz) , 
the  period  used  for  time-domain  mb,  the  number  of  the  zero 
crossing  used  for  the  mb  measurement  (b  indicates  measurement  at 
270  degrees  azimuth),  and,  finally,  the  time-domain  mb 
measurement.  In  general,  the  time-domain  mb  measurements  are 
larger  than  the  frequency-domain  estimates.  There  are  three  reasons 
for  this.  First,  the  t*  of  .8  may  be  lower  than  t*  of  NTS.  Bache 
et  al.  (1980)  estimate  a  t*  of  .8  to  1.05  for  NTS  explosions 
observed  at  different  receiver  stations.  Hadley  (1979)  estimates  t* 
for  NTS  to  lie  in  the  range  of  1.1  to  1.5  (Both  t*  estimates  were 
for  a  frequency  of  about  one  Hertz).  Since  our  frequency-domain 
estimates  are  caliDrated  based  on  NTS  shots  in  tuff,  this  increases 
both  explosion  and  earthquake  magnitudes  by  about  a  quarter  of  a 
magnitude  compared  to  using  a  t*  of  1.0.  Another  possibility  is 
tnat  the  mantle  geometric  spreading  factor  used  here  is  too  large. 
Second,  the  orientation  of  the  faults  puts  the  earthquake  rays 
closer  to  a  radiation  pattern  maximum  than  the  average  used  in 
making  frequency-domain  estimates.  This  increases  the  amplitudes  of 
the  earthquakes  by  another  quarter  magnitude.  Finally,  there  is  an 
increase  in  amplitude  caused  by  the  shift  in  the  dominant  frequency 
away  from  the  peak  when  the  peak  instrument  response  is  above  the 
corner  frequency.  Notice  tne  increase  in  period  as  a  function  of 
the  moment  of  the  earthquake. 

Figure  17  is  a  plot  of  time-domain  mfa  measurements  for  all 
of  tne  earthquake  and  explosion  sources  listed  in  the  tables. 
Explosion  synthetics  were  made  at  1  KT,  10  KT,  100  KT  and  1000  KT. 
This  figure  may  be  compared  with  the  frequency-domain  figure  (Figure 
9)  comparing  mb  and  M$.  The  separation  between  the  earthquake 
and  explosion  populations  is  still  very  clear;  however,  they  appear 
to  converge  at  small  magnitudes.  The  reason  for  this  is 
interference  with  pP.  For  the  shallow  1  KT  event  in  granite, 
interference  reduces  mb  by  about  0.7;  at  10  KT  it  reduces  mb  by 
0.2;  at  100  KT  it  has  little  effect  while  at  1000  KT  it  increases 
mb  by  about  0.3. 
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Figure  17.  Time-domain  mb  measurements  from  synthetic  seismograms 
of  the  six  earthquake  sources  scaled  to  different  sizes 
and  Mueller-Murphy  explosion  sources  in  salt  (S),  tuff 
(T),  granite  (G),  and  shale  (H)  at  1,  10,  100  and  1000 
KT.,  Convergence  at  low  magnitude  is  caused  by  increas¬ 
ing  destructive  interference  of  P  and  pP  as  explosions 
decrease  in  depth  with  decreasing  yield. 
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VIII.  RECOMMENDATIONS  FOR  FURTHER  RESEARCH 


In  this  report,  we  have  tried  to  identify  and  quantify  the 
important  effects  which  allow  or  hinder  discrimination.  It  is 
impossible,  however,  to  examine  every  characteristic  of  earthquakes 
or  explosions  which  might  apply  to  this  problem.  We  would  like  to 
mention  some  of  the  factors  which  we  have  not  discussed  in  detail. 

The  extent  of  pP  interference  is  quite  important  to  both  VFM 
and  m|jiMs  discrimination,  especially  for  small  explosions. 
Further  work  on  how  near  surface  nonlinearities  such  as  cracking  and 
spallation  affect  pP  is  needed.  This  is  an  active  area  of  research, 
and  nonlinear  effects  of  pP  are  currently  being  modeled  by  two 
dimensional  (cyl indrical ly  symmetric)  nonlinear  numerical 
calculations  (Oay,  et  jT. ,  1981). 

Our  analysis  of  the  frequency  distribution  of  m^  and  M$ 
was  performed  for  a  uniform  half-space  and  then  extended  to  a  layer 
over  a  half-space.  It  would  be  possible  to  perform  similar 
calculations  for  full  layered  half-space  models.  There  would  then 
be  variations  in  seismic  radiation  due  to  the  depth  dependence  of 
source  region  elastic  properties  and  differences  in  the 
eigenfunctions  of  the  layered  half-space. 

Von  Seggern  and  Rivers  (1979)  performed  an  analysis  similar  to 
ours  using  a  Gutenberg  earth  model.  Their  use  of  the  root  mean 
square  amplitude  of  the  distribution  function  rather  than  the  peak 
of  the  distribution  function  (Figure  4)  decreases  the  effect  of 
source  mechanism  on  earthquake/explosion  separation  from  .35  to  .22 
magnitude  units.  They  also  find  this  number  to  be  consistent  in  the 
upper  40  km  except  for  a  decrease  to  zero  over  a  narrow  range  of 
depths  near  5  km. 

The  best  approach  to  the  layered  space  problem  would  oe  to 
evaluate  the  relative  excitation  of  surface  waves  for  individual 
structures  of  interest.  It  would  also  be  possible  to  perform  the 
analysis  for  a  restricted  range  of  depths  and/or  orientations  if  the 
earthquakes  in  a  given  region  are  known  to  have  preferred  depths  or 
orientations. 


The  six  earthquake  models  in  Figure  1  are  based  on  plausible 
physical  assumptions,  and  exhibit  rupture  histories  and  average 
stress  drops  which  are  consistent  with  what  is  known  observational ly 
aoout  natural  earthquakes.  Nonetheless,  they  may  not  be  fully 

representative  of  the  entire  range  of  phenomena  relevant  to  the 
discrimination  problem.  One  case  of  potential  interest  not 
considered  here  is  that  of  a  fault  with  a  localized  stress 
concentration  many  times  greater  than  the  average  stress  drop.  Such 
a  configuration  has  been  inferred  for  the  1971  San  Fernando 
Earthquake  by  several  authors  (e.g.,  Bache  and  Barker,  1978;  Heaton 
and  Helmberger,  1979).  Some  insight  into  this  question  is  provided 
by  Models  1  and  2.  We  find  that  Model  2  is  easier  to  discriminate 
than  Model  1  in  spite  of  the  fact  that  both  have  the  same  stress 
concentration  at  the  center,  indicating  that  an  earthquake 
containing  a  stress  concentration,  but  which  propagates  into  a 
region  of  lower  stress,  will  be  easier  to  discriminate  than  failure 
of  an  isolated  stress  concentration  which  does  not  propagate  further. 

Our  analysis  of  the  VFM  discriminant  was  based  on  spectral 
amplitudes.  In  practice,  measurements  are  made  in  the  time  domain 
after  narrow  band  filtering,  and  restricted  to  the  first  observed 
arrival  at  each  frequency.  There  are  also  several  different  ways  of 
smoothing  the  spectrum.  It  would  be  valuable  to  apply  current  VFM 
methods  to  a  range  of  synthetic  seismograms  from  earthquake  and 
explosion  models. 

Another  effect  which  can  be  important  is  explosion  induced 
tectonic  strain  release.  The  effect  of  tectonic  release  on  body  and 
surface  waves*  is  discussed  in  detail  using  a  linear  approximation 
for  explosions  In  complex  prestress  fields  by  Stevens  (1982).  In 
general,  tectonic  release  adds  energy  to  the  seismic  radiation  and 
imposes  a  radiation  pattern  on  the  body  and  surface  waves.  Tectonic 
release  can  make  discrimination  more  difficult,  particularly  if 
observations  cover  only  a  limited  portion  of  the  focal  sphere.  In 
some  cases,  tectonic  release  can  aid  discrimination  by  reducing 
surface  wave  amplitudes  and  increasing  high  frequency  body  wave 
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amplitudes.  Nonlinear  simulations  of  explosion  induced  tectonic 
release  (e.g.,  Oay,  et  _a1_. ,  1981)  are  Important  in  order  to  better 
determine  the  tectonic  source  and  to  understand  how  it  interferes 
with  and  influences  the  explosion  source. 

Finally,  we  note  that  all  of  the  analysis  was  performed  for 
earthquakes  and  explosions  in  similar  structures  and  may  therefore 
not  apply  to  the  comparison  of  earthquakes  in  one  area  with 
explosions  in  another.  Variations  in  attenuation  and  the  frequency 
dependence  of  attenuation,  as  well  as  the  amplitude  response  of  a 
given  structure,  may  strongly  affect  the  discriminants.  These 
variations  can  only  be  determined  by  a  careful  study  of  structure 
and  attenuation  in  all  areas  of  interest. 
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IX.  DISCUSSION  ANO  CONCLUSIONS 

The  theoretical  results  predict  that  both  m^:.^  and 
variable  frequency  magnitude  methods  will  successfully  discriminate 
most  earthquakes  from  underground  explosions,  but  will  be  affected 
in  very  different  ways  by  a  number  of  important  factors.  While 
source  spectral  differences  are  important  to  both  discriminants, 
differences  in  source  mechanism  (center  of  dilatation  versus 
double-couple)  and  differences  in  source  region  elastic  properties 
may  dominate  for  small  events  using  the  m^:Ms  discriminant. 
Free-surface  reflected  phases  also  play  an  important  role  at  low 
magnitude,  tending  to  innibit  discrimination  and  enhance 

VFM  discrimination.  Interference  by  pP  is  sufficient  to  cause  the 
theoretical  explosion  population  to  begin  to  converge  with  the 
earthquaxe  population  on  the  .  mb:M$  plane  at  about  mb  *  4.5. 
The  latter  conclusion,  of  course,  depends  on  the  assumption  that  the 
explosions  are  burieo  at  normal  containment  depths  and  that  pP  is 
approximately  elastic. 

The  average  separation  between  the  theoretical  earthquake  and 
explosion  populations  on  the  mb:M$  plane,  excluding  interference 
effects  due  to  surface  reflected  phases,  is  given  approximately  by 


.35  ♦  .5 


+  spectral 
differences 


» 


where  x  indicates  explosion  and  q  indicates  earthquake,  and  o  and  a 
are  source-region  density  and  P-wave  speed.  The  first  term 
represents  tne  effect  of  source  mechanism,  which  Is  responsible  for 
an  average  m^Mj  separation  of  .35  magnitude  units  between 
earthquake  and  explosion  populations.  For  an  earthquake  and 
explosion  of  equal  seismic  moment,  mb  of  the  explosion  exceeds  the 
average  earthquake  mb  by  .45  units,  while  M$  of  the  explosion 
exceeds  the  average  earthquake  Ms  by  .1  unit.  As  shown  in  Figures 
4  and  5,  considerable  scatter  is  Introduced  by  variations  in  the 
source  mechanism. 
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The  second  term  represents  the  effect  of  source-region  elastic 
properties.  Source-region  properties  will  usually  increase  the 
separation  between  earthquakes  and  explosions  on  the  mb:Ms  plane 
since  explosions  are  shallower  than  most  earthquakes,  and  so  tend  to 
be  in  lower- impedance  materials.  Small  explosions  in  low-velocity, 
low  density  material  are  predicted  to  discriminate  much  more  easily 

oy  mb:Ms  than  are  small  explosions  in  high-velocity, 

high-density  material.  This  is  true  even  though  the  explosion 
source  spectrum  is  generally  more  earthquake-like  (lower  comer 
frequency,  less  peaked)  in  soft  materials  than  in  hard  materials  as 
can  be  seen  in  Figure  11,  for  example.  In  other  words,  in  our 
theoretical  results,  the  elastic  wave-propagation  effect  of  the 
source  medium  dominates  over  its  competing  effect  on  the  explosion 
source  spectrum  for  m^  less  than  about  5.5  (at  one  Hertz). 

The  third  term,  source  spectral  differences,  accounts  for  most 
of  the  increasing  separation  predicted  above  mb=s  4.5 

(see  Figures  9  and  17).  A  consistent  feature  of  the  spectral 
magnitude  estimates  (Figure  8)  and  the  synthetic  seismogram 
magnitude  estimates  (Figure  16)  is  that  Earthquake  Model  3  falls 
nearest  the  explosion  population  followed  closely  by  Models  4  and  1, 
while  Models  2,  5,  and  6  fall  farther  away.  A  crude  explanation  may 
oe  given  for  this  behavior  by  considering  three  factors:  rupture 
symmetry,  rupture  duration,  and  stress  drop.  The  difficulty  of 
discrimination  is  given  approximately  for  each  event  by  the 

3 

(constant)  product  of  MQfc  where  Mq  is  the  moment  and  fc 

is  the  corner  frequency.  The  corner  frequency  is  (approximately) 
inversely  proportional  to  the  rupture  duration;  so  a  slow  rupture 
velocity  such  as  in  Model  2  tends  to  enhance  discrimination.  In 
asymmetric  faults  such  as  Models  5  and  6,  the  fault  area,  and 
therefore  the  moment,  is  reduced  when  compared  to  a  symmetric  model 
such  as  Model  3,  but  the  rupture  duration  is  still  determined  by  the 
longest  dimension  so  that  the  corner  frequency  is  approximately 
unchanged.  Again,  this  enhances  discrimination.  Finally,  the 
spontaneous  rupture  models  (Models  1,  2,  and  6)  all  extendea  into 
regions  of  lower  prestress  so  that  the  average  stress  drops,  and 
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therefore  the  moments,  of  these  events  are  somewhat  lower  than  the 
other  mooels.  This  is,  of  course,  a  very  rough  explanation  which 

*  combines  subtle  details  of  the  amplitude  and  phase  spectra  of  the 
source  into  the  single  gross  feature  "corner  frequency,"  but  it 
provides  a  partial  physical  explanation  of  the  variability  in 
discrimination  potential  between  the  earthquake  models. 

The  important  result  here  for  di scrimination  studies  is  that 
the  most  difficult  earthquakes  to  discriminate  are  those  which  are 
symmetric  (approximately  square  or  circular),  have  rapid  rupture 
growth,  and  have  uniformly  high  stress  drop.  This  will  be  true  both 
for  mj?:Ms  and  VFM  discrimination.  An  interesting  fact  is  that 
the  complexity  of  Model  6  significantly  enhances  the  predicted 
discrimination.  Also,  note  that  Model  2  is  easier  to  discriminate 
than  Model  1  in  spite  of  the  fact  that  they  both  have  the  same 
stress  concentration  at  the  center. 

Variable  frequency  magnitude  discrimination  depends  on 
spectral  differences  between  earthquakes  and  explosions.  The  same 
I  spectral  features  of  the  earthquake  models  which  affect  mb:Ms 

also  affect  VFM.  Models  2,  5,  and  6  separate  from  explosions  more 
easily  than  Models  3,  4,  and  1  do.  In  most  cases,  we  find  that  the 
theoretical  explosion  source  function  has  more  high-frequency  energy 
j  than  tne  eartnquaxe  source  functions  (scaled  to  100  bars  maximum 

stress  drop);  so  the  theoretical  event  populations  are  separated 
down  to  about  «  3.5.  Explosions  in  saturated  granite  separate 
at  even  smaller  magnitudes. 

*  Explosions  in  high  velocity  materials  such  as  granite  should 

be  more  readily  identified  by  VFM  than  explosions  in  lower  velocity 
materials  such  as  tuff,  which  is  opposite  to  the  low-magnitude 
behavior  of  mb:Ms<  For  small,  shallow  explosions  in 

I  low-velocity  materials  such  as  tuff,  source  spectral  differences 

alone  may  be  insufficient  to  achieve  separation  from  the  earthquake 

population  (see  Figure  13).  Interference  of  pP  can  lower  the 
low-frequency  magnitude  of  shallow  explosions  by  up  to  half  a 
»  magnitude  unit  and  should  be  adequate  to  permit  most  such  events  to 
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be  discriminated.  However,  we  have  not  performed  a  synthetic 
seismogram  study  of  the  VFM  discriminant  to  confirm  this 

*  conjecture.  Conversely,  over-burial  of  small  explosions  in  soft 
rode  may  defeat  the  VFM  discriminant.  Of  course,  deep  burial  in 
soft  rock  should  favor  event  identification  by  mb:Mg. 

To  a  surprising  degree,  we  have  found  that  the  mb:M$  and 
VFM  methods  complement  each  other.  The  VFM  method  is  predicted  to 
be  most  effective  for  mb  ranging  from  about  3.5  to  5.5  whereas 
mb:M$  is  predicted  to  be  most  effective  above  about  m^  «  4.5 
and  to  degenerate  as  magnitude  decreases  below  mb«4.5,  due  to  pP 
interference.  Furthermore,  at  low  magnitudes,  mb:Ms  is  favored 
by  soft  explosion  media,  and  VFM  is  favored  by  hard  explosion 
media.  As  a  third  example,  we  find  discrimination  of  small  events 

by  m.:M  to  be  virtually  uninfluenced  by  earthquake  stress-drop 

i  D  S 

whereas  the  VFM  method,  depending  as  it  does  primarily  on  source 
spectral  difference,  could  conceivably  be  defeated  by  earthquakes 
with  anomalously  high  stress-drops.  Finally,  mb:Ms  fails  for  a 
few  earthquake  source  orientations  and  depths  for  which  surface  wave 
excitation  is  anomalously  low  (notably  very  shallow,  vertical 
dip-slip  faults)  while  VFM  is  only  slightly  affected  by  source 
orientation.  On  the  basis  of  this  study,  therefore,  we  conclude 
that  simultaneous  application  of  the  mb:Ms  and  VFM  methods  will 

*  substantially  increase  the  likelihood  of  correct  event 
identification. 
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APPENDIX  A 

MUELLER-MURPHY  EXPLOSION  SOURCE  FUNCTION 
AND  SCALING  LAWS 


The  Mueller-Murphy  empirical  scaling  laws  give  the  source  time 
function  for  explosions  In  different  materials  as  a  function  of 
yield  and  depth  (Mueller  and  Murphy,  1971;  Murphy  and  Mueller,  1971; 
Murphy  1977a,  1977b).  As  a  general  rule,  larger  explosions  will  be 
more  deeply  buried  (for  containment).  We  expect  depth  to  scale 
approximately  with  the  cube  root  of  yield;  that  is. 


> 


where  h  is  the  explosion  depth,  W  is  the  yield  and  h  is  an 
empirical  constant.  For  most  materials  hQ  »=  122  m/KT1  .  The 
source  is  taken  to  be  a  pressure  pulse  of  the  form 

P(t)  -  {p0e'7t  +  V  H{t)  ’ 

applied  at  the  elastic  radius  R.  The  values  of  PQ,  P2,  y  and  R 
are  given  by 


P1  =  P1Q  h/h0 

R  *  R„  (h./h)1/n  W1/3 

o  o 

P2  =  P2Q  (Vh)1/3  (VR)3  W’8? 

t  *  rQ  R0/R 

Po  aPl-P2 

The  empirical  constants  are  given  in  Table  A  for  yield  in  kilotons 
(Murphy,  1977b).  The  reduced  velocity  potential  (i^)  is  then  given 
by: 
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Hu) 


e1uR/a 


u  +  1  u  u  - 
0  0 


X  +  2  u 


2 

u 


where  a  is  the  P-wave  velocity,  x  and  u  are  the  elastic  constants 
and  uQ  *  a/R. 

The  far-field  displacment  (for  any  RVP)  is  given  at  a  distance 

r  by 


u(r,u) 


•  jur 

'P(u)  e  a 

ar 


The  long  period  limit  'P „  is  given  by 


P20  Ro  ,,.76 

3U  w 


for 


hoM 


1/3 


1/3 

For  all  materials  except  salt,  if  h  *  trW  ,  then  P,  scales 

1/3  0  L 

as  W  while  P7  and  R,  which  are  affected  by  overburden 

4  20 
pressure,  vary  approximately  as  W'fc  . 
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MUELLER/MURPHY  Et-PIRICAL  CONSTANTS 


Constants 

Tuff/Rhyolite 

Granite 

Salt 

Shale 

ho 

122 

122 

122 

122 

meters 

Ro 

202 

321 

478 

265 

meters 

To 

26 

34 

31 

42 

sec"1 

P10 

36 

46 

55 

42 

bars 

P20 

50 

24 

8 

25 

bars 

n 

2.4 

2.4 

1.87 

2.4 
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APPENDIX  B 


EQUIVALENT  JUMP  VECTORS  FROM  UPGOING  AND  DOWNGOING  WAVES 


We  want  to  make  synthetic  seismograms  using  the  method  of 
Bache  and  Harkrlder  (1976)  when  the  source  Is  described  In  terms  of 
upgolng  and  downgoing  waves.  Bache  and  Harkrlder's  method  starts 
with  a  multipolar  expansion  for  the  source,  derives  equivalent  jump 
vectors  and  uses  a  saddle  point  approximation  for  the  far  field 
displacement. 

In  general,  a  source  may  be  described  in  cylindrical 
coordinates  by  three  potentials: 

P(r,  e,  2,  u)  *  I  f  si  (u,  k)  e‘i,cra'z"h I  J  (kr)  dk 
mt0  JQ  r,n  m 

X'ir,  e,  z,  u)  *  £  f  Si  (u,  k)  e~1krV2'h*  J  (kr)  dk 
m»o  ■'o 

oo  oo 

X  ( r,  e,  z,  u)  -  I  s^tu,  k)  e'ikr8l2'hl  Jm(kr)  dk 

where  +  correspond  to  z  %  h,  with  h  the  source  depth.  The 
potential  X  corresponds  to  SH  waves  and  is  decoupled  from  the  other 
potentials.  Using  the  saddle  point  approximation  as  in  the  Bache 
and  Harkrider  paper,  the  far-field  displacements  are  found  to  be 


U 

U 

U 


P 

SV 

SH 


.  _  -ik  R 

I  St  1mr  k 
~  itn  a  a  K 


_  y  e+  ^m+l  _  1,1,  e _ s 


>  i .  r„kk 
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-Ik.R 
s**b  IT 


ik.R 


-2 
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we  have  used  the  definition 


r 


1/2 

1) 


where  c  is  the  phase  velocity-  In  the  far-field  approximation, 

y 

k  =  u/c  corresponds  to  a  particular  ray  for  which  c  -  --j-—  where  e 
is  the  takeoff  angle  of  the  ray  and  v  is  the  appropriate  velocity  a 
or  s-  The  equation  for  (in  Equation  33  of  Bache  and 

Harkrlder)  is  in  error  and  should  be  multiplied  by  -1. 

We  compare  these  equations  with  the  far-field  displacements 
for  the  earthquake  source: 


V 


usv- 


'SH 


1 

-Ik  R 

e  a 

7  3“ 

4irpa 

ft 

1 

-Ik  R 
e  8 

4irp83 

ft 
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-ik.R 
e  3 

4»pS3 
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Rp(e,  4)  ftp(e,  4, 
4)  ^^(e*  d» 


Sj^(e,  4)  ns(e,  4, 
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U 


u) 


) 

) 


Here  R  ,  Rc  ,  Rc,  are  radiation  pattern  coefficients  (see 

P 

Appendix  0)  which  are  independent  of  frequency,  and  and 
are  source  spectra  obtained  by  integrating  the  slip  velocity  over 
the  fault  plane  (see  Appendix  C).  Comparing  these  two  sets  of 
equations,  we  find 


Rq(«,  4)  ap  (e,  4,  u) 

4*pa  r  k 
a  a 

1R$^(0f  tf)  u) 

4,p8  krBkB 

-Rj^  (®,  4)  d, 

4ir0e3  krg 
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where  e,  i  and  k  all  correspond  to  the  appropriate  directions  for 
upgoing  and  downgoing  waves. 

From  this  result  we  can  derive  equivalent  jump  vectors  for  the 
earthquake  source.  Using  the  notation  of  Bache  and  Harkrlder  and 
the  results  of  Harkrider  (1964), 


sUm  *  k2[(Slm  *  sIm'  -  ,kVs+2m  * 

•».  *  ikt-""-.'sI,,+  SI»' +  k2<4  - 

.z.-kVpCIt-iks^-  sj.)  ,1] 

-  k2c2oC-Tra(S*„,  *  Sjm)-ik  (,  -  1)(S*„  -  SJJ] 


where 

.  2 

y  ■  2  <t>  • 

For  the  tk-gential  waves 


* ik2  (SL  - 


-  -,k2“r,<S*3™  *  S‘s,» 


The  equivalent  jump  vectors  in  terms  of  upgoing  and  downgoing  waves 
are  defined  by 

iV  «  <sU_ 


with  similar  definitions  for  the  other  jump  vectors.  These  may  be 
evaluated  directly  from  the  upgoing  waves  and  used  in  Equation  28  of 
Bache  and  Harkrider.  This  removes  the  m  index  on  all  remaining 
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1 


equations  with  the  final  result  for  the  far-field  displacement  from 
an  earthquake  In  a  layered  medium  {equivalent  to  their  Equation  33) 
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APPENDIX  C 


INTERPOLATION  METHOD 


The  evaluation  of  the  source  time  function  from  the  finite 
difference  calculations  requires  the  evaluation  of  surface  integrals 
over  the  failure  surface  which  can  be  reduced  to  a  three-dimensional 
Fourier  transform  (Day  et  al . ,  1978).  The  integrals  have  the  form 


ap(e,rf,u)  ■  u  J  J  J  dxQdy0dt  i(xQ,yQ,t)  ei(Kxx*V”ut) 


fls(e,rf,u)  ■  u  J  f  /  dxQdyQdt  i(x0,yQ,t)  el(lV»+Kyy"ut) 


K  *  —  sine  cosd 

X  a 


Ky  *  «  s-ine  sin** 


K  *  sine  cosd 
x  b 

Ky  *  j  sine  sind 


where  s  is  the  slip  velocity  which  is  known  at  evenly  spaced 
Intervals  of  xQ,  yQ,  and  t.  The  x  and  y  coordinates  are 

Cartesian  coordinates  in  the  fault  plane.  The  angles  e  and  i  are 

spherical  coordinates  of  the  observer,  with  e  measured  from  the 

fault  normal,  and  i  measured  in  the  fault  plane. 

The  results  of  the  surface  integration  give  a  ,  at 

evenly  spaced  intervals  of  K  ,  K  when  an  FFT  is  used.  For  the 

a  y 

computation  of  a  particular  ray  leaving  at  an  angle  e,  «f,  we  will 
therefore  need  to  Interpolate  in  the  K  ,  K  plane  to  find  the 

*  J 

spectrum  at  the  appropriate  ray  parameters  for  each  frequency. 

The  method  of  Interpolation  is  described  by  Bracewell  (1965) 
and  consists  of  convolution  of  the  function  known  at  evenly  spaced 
points  with  the  two-dimensional  sine  function 


S(IW 


*Ky 

sin-^ 

_ y 
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The  interpolated  values  will  be  exact  as  long  as  the 
transformed  function  is  zero  outside  of  the  interval 


<  x  < 


-V*<T* 


where  x,  y  are  the  transform  variables  correspondi ng  to  K  ,  K  . 

a  y 

We  can  write  the  function  at  discrete  points  in  the  form 

OO  oo 

f I nt  .  m  „ 

y 


F{Kx,  K  )  *  I  I  mV  s(Kx"nTx)  5(Ky"mV 

■>  n=-oom=-oo  J  j  j 


Convolution  with  the  sine  function  gives 


f(K  ,K  )  *-£*?•[[ dx*  dy ’  I  f(nx  ,mx  )  s(x'-nx  )  «(y'-nx  ) 
x  y  JJ  m,n  *  ■>  *  •> 


si  rife-  (K  -x*)]  sinfe  (K  -y')] 
- (*.-»'>  I  K  -/)  ^  ^ - 


If  Kx  a  qxx  and  »  pTy,  then,  using 


sin(irq  -  mr)  *  (-l)n  sinirq  , 
we  have 


*/„  lr  *  _  slrnrq  sirnrp  v  ,  nn+m  f(nTx,mV 

f(Kx,Ky) - y -  l  (-1)  Iq-nKp-ii 


If  either  p  or  q  is  an  integer,  the  integral  becomes 
one-dimensional.  If  Kx  *  nxx,  then 

,{n.n  r  m  f(nxv,mx  ) 
f(Kx.Ky)  .  UffiE  2  l-l)m  X 


m 


p-m 


The  one-dimensional  interpolator  is  also  very  useful  for 
interpolating  the  frequency  spectrum.  If  u  *  px,  then 


f(u)  -  S  (-1)"  liail 


81 


scasEO 


1 


This,  again,  will  be  exact  as  long  as  the  transformed  function 
exists  only  in  the  time  interval 


1 

'  It 


<  t  < 


1 

It 
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APPENDIX  D 

UPGOING  AND  DOWNGOING  'WAVES  FOR  ARBITRARY  FAULT  ORIENTATION 
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APPENDIX  D 


UPGOING  AND  DOWNGOING  WAVES  FOR  ARBITRARY  FAULT  ORIENTATION 

As  described  in  the  previous  appendix,  a  surface  Integral  can 

be  performed  to  give  the  outgoing  waves  in  any  direction.  We  are 

Interested  in  the  waves  which  correspond  to  upgoing  and  downgoing 

waves  in  a  half-space.  The  source  and  half-space  coordinate  systems 

are  shown  In  Figure  18.  We  need  to  determine  the  angles  e  and  d  in 

the  fault  coordinate  system  for  both  the  upgoing  (-)  and  downgoing 

(+)  waves  for  an  arbitrary  fault  orientation.  These  angles  will 

then  be  used  to  find  the  appropriate  values  of  K  and  K  as 

a  y 

described  in  the  last  appendix  and  to  include  the  radiation  pattern. 

From  Aid  and  Richards,  the  slip  and  normal  vectors  are  given 
in  terms  of  the  strike  d$,  dip  5,  and  rake  x  in  the  half  space 
coordinate  system  by 

s  a  (cosx  cosd  +  cosa  sinx  sind  » 
s  s 


cosx  slnd  -  cosa  sinx  cosd  ,  -sinx  slnr) 
s  s 


(D.l) 


n  *  (-  sina  sind  ,  sina  cosd  ,  -  cosa) 

s  s 


The  vectors  K+  and  K"  are  defined  by  the  takeoff  angle  and 
azimuth  of  the  observation  point 

K+  *  (sinet  cosdt,  sinet  sind^,  coset) 

K"  -  (sinefc  cosdt,  sine*  sindt>  -  coset)  (D.2) 
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Figure  18 


Fault 

Coordinate 

System 
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Fault  and  global  coordinates  used  for  synthetic  seismograms. 
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From  these  we  need  the  angles  In  the  earthquake  frame  e~, 

A  unit  vector  In  the  earthquake  coordinate  system  is  given  by 

e  *  (sine  cosrf,  sine  slnrf,  cose) 

with 

s  *  (0,  1,  0),  n  *  (0,  0,  1)  and  s  x  n  3  (1,  0,  0)  ; 

A 

so,  using  K  as  the  unit  vector, 
n  X-  *  +  cose- 
s  *  X—  »  sine^  sini£ 

(s  x  n)  *  X—  3  sinei  cos^- 


Using  Equations  0.1  and  0.2  to  compute  n,  s,  X+  and  X"  in  the 

half  space  frame  gives  the  angles  e*,  t*  in  the  earthquake 

frame.  The  values  of  X  ,  X  can  now  be  found: 

x  y 

For  P-waves 

xi  3  -  sine^-  cos^ 

X  21 

X*  =*  -  sinei  s  in«£ 

y  a 

For  S-waves 

t  sine^  costf^ 
x  a 

X j  j  sine^  sin^ 

For  a  particular  problem,  the  same  X*,  X*  is  used  for  both 

x  y 

P-waves  and  S-waves.  e^ ,  ^  correspond  to  the  takeoff  angle  of  the 
wave  of  interest. 


I 


» 
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The  radiation  pattern  factors  in  the  earthquake  coordinate 
system  are 


sin2s  si nd 
cos2e  sind 
cose  cosii 
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